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Detachable  Summary 

The  Federal  Emergency  Management  Agency  (FEMA)  has  supported  a  contract  in  which  the  Law¬ 
rence  Livermore  National  Laboratory  (LLNL)  has  acted  as  contractor  and  Science  Applications,  Inc.  (SA1) 
the  subcontractor  to  develop  a  coupled  set  of  models  between  those  of  the  LLNL  Atmospheric  Release 
Advisory  Capability  (ARAC)  system  and  candidate  evacuation  models. 

A  rapid  computer  code  easily  adapted  for  use  with  color  graphics  was  developed  to  convey  the 
dynamics  of  population  evacuation  in  the  vicinity  of  a  nuclear  power  generating  station.  The  model 
developed  required  the  evacuation  model  to  integrate  calculated  meteorological  transport  and  diffusion  of 
an  atmospheric  release  and  the  resulting  dose  commitments  as  a  function  of  time  and  position  from  the 
release  point.  An  integral  part  of  the  demonstration  is  the  display  of  isopleths  resulting  from  an  airborne 
radiactive  release  using  ARAC  data  to  show  potential  evacuation  problems. 

The  ARAC  system  utilizes  a  suite  of  numerical  models  appropriate  to  a  variety  of  atmospheric  release 
incidents.  For  the  purpose  of  this  study,  concentration  contours  coupled  with  the  SAI  evacuation  model 
were  calculated  by  using  the  MATHEW  and  ADPIC  codes. 

The  evacuation  emergency  response  model  coupling  ARAC  output  was  developed  for  use  on  a  VAX- 
780  in  conjunction  with  a  VSV11  color  graphics  terminal.  LLNL  designed  the  graphics  package.  Informa¬ 
tion  is  available  to  users  of  the  model  as  a  function  of  time  to  assist  in  either  training  sessions  or  an  actual 
emergency  response.  At  the  beginning  of  an  accident  and  for  every  15  minutes  of  simulated  time  there¬ 
after,  the  user  may  do  the  following: 

•  View  a  colored  map  indicating  the  road  network  and  the  number  of  individuals  within  a  10-  and 
15-mi  radius  of  the  accident  at  the  NPP. 

•  Display  dose  contours  for  1311, 133Xe,  and  137Cs  overlayed  on  the  colored  map. 

•  Display  a  histogram  showing  the  number  of  people  still  within  the  evacuation  zone,  those  evacu¬ 
ated,  and  accumulated  population  doses. 

•  Determine  directly  from  the  color  graphics  displayed  whether  a  bottleneck  exists  in  evacuating 
individuals  along  a  particular  clogged  road. 

•  Determine  directly  from  the  color  graphics  and  overlayed  contours  whether  a  particular  evacua¬ 
tion  routefs)  intersects  the  radiation  plume  endangering  individuals  being  evacuated. 

•  Skip  any  of  the  preceding  options  or  backtrack  to  further  study  a  particular  display. 

The  EVACD  computer  program  was  developed  to  demonstrate  the  feasibility  of  coupling  an  evacua¬ 
tion  model  with  a  color  graphics  system  to  illustrate  the  useful  information  available  from  such  a  com¬ 
bination.  The  objective  of  the  program  was  achieved  and  has  enabled  improvement  of  the  system.  The 
following  are  suggestions  for  improvement  to  the  developed  model: 

1.  Add  more  interactive  capabilities  to  allow  operators  running  the  model  to  modify  road  network 
conditions  at  each  appropriate  time  step.  In  other  words,  the  operator  should  be  allowed  to  insert 
road  blocks,  change  road  preference  factors,  and  similar  types  of  real-time  changes. 

2.  Add  more  graphics  to  make  the  model  the  basic  tool  for  training  emergency  procedure  personnel. 
Suggested  output  would  aid  in  assessing  the  various  evacuation  decisions. 

3.  Generalize  the  dose  assessment  routines  to  cover  more  cases  of  a  general  nature. 

4.  Refine  the  road  network  model  to  include  more  specific  data,  such  as  road-dependent  capacities. 

5.  Calibrate  the  model  against  data  for  past  actual  evacuations. 

6.  Develop  a  training  program  using  EVACD  as  the  main  teaching  tool. 
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Abstract 

Federal  Emergency  Management  Agency  (FEMA)  sponsored  project  to  develop  a 
coupled  set  of  models  between  those  of  the  Lawrence  Livermore  National  Laboratory 
(LLNL)  Atmospheric  Release  Advisory  Capability  (ARAC)  system  and  candidate  evacua¬ 
tion  models  is  discussed  herein.  LLNL  and  Science  Applications,  Inc.  (SAI)  serve  as  con¬ 
tractor  and  subcontractor,  respectively. 

This  report  describes  the  ARAC  system  and  discusses  the  rapid  computer  code  devel¬ 
oped  and  the  coupling  with  ARAC  output.  The  computer  code  is  adapted  to  the  use  of 
color  graphics  as  a  means  to  display  and  convey  the  dynamics  of  an  emergency  evacua¬ 
tion.  The  model  is  applied  to  a  specific  case  of  an  emergency  evacuation  of  individuals 
surrounding  the  Rancho  Seco  Nuclear  Power  Plant,  located  approximately  25  miles  south¬ 
east  of  Sacramento,  California.  The  graphics  available  to  the  model  user  for  the  Rancho 
Seco  test  case  are  displayed  and  noted  in  detail.  Suggestions  for  future,  potential  improve¬ 
ments  to  the  emergency  evacuation  model  are  presented. 


Introduction 


The  Federal  Emergency  Management  Agency 
(FEMA)  has  supported  a  contract  in  which  the 
Lawrence  Livermore  National  Laboratory  (LLNL) 
has  acted  as  contractor  and  Science  Applications, 
Inc.  (SAI)  the  subcontractor  to  develop  a  coupled 
set  of  models  between  those  of  the  LLNL  Atmo¬ 
spheric  Release  Advisory  Capability  (ARAC)  sys¬ 
tem  and  candidate  evacuation  models.  The  LLNL 
ARAC  staff  has  worked  in  conjunction  with  the 
staff  of  SAI  to  meet  the  following  objectives: 

1.  Couple  the  atmospheric  and  dose  models 
used  in  the  ARAC  system  with  evacuation 
models  developed  by  SAI  as  part  of  a  study 
for  the  California  Office  of  Emergency 
Services. 

2.  Study  the  feasibility  and  effectiveness  of 
putting  the  output  of  evacuation/ 


emergency  response  models  such  as  the  lat¬ 
est  U.S.  Department  of  Transportation 
package  UTPS  (Urban  Transportation  Plan¬ 
ning  System)  or  other  network  models  of 
choice  into  color  graphic  output  format 
compatible  with  ARAC's  CRT  outputs. 

3.  Integrate  evacuation  models  into  the  ARAC 
services  package,  thereby  extending  the  ca¬ 
pability  of  the  ARAC  system. 

4.  Test  and  evaluate  the  concept  of  integrating 
dose  and  evacuation  models  utilizing  data 
developed  for  the  Rancho  Seco  Nuclear 
Power  Station  site  in  the  State  of  California 
by  means  of  a  working  demonstration  using 
available  equipment. 


The  Atmospheric  Release  Advisory  Capability  (ARAC) 


Introduction 

The  ARAC  project  was  initiated  by  the 
Atomic  Energy  Commission,  now  the  Department 
of  Energy  (DOE),  to  develop  a  computer-based 


capability  to  produce  rapid  projections 
(advisories)  of  the  transport,  diffusion,  and  depo¬ 
sition  of  radioactive  material  released  into  the  at¬ 
mosphere.  The  concept's  feasibility  was  demon¬ 
strated  in  1975,  after  which  time  the  central 


facility  was  established  at  LLNL  and  computer 
data  processing  equipment  was  installed  at  se¬ 
lected  DOE  sites.1 

The  ARAC  system  (Fig.  1)  makes  available  to 
users  predictive  data  from  proven  and  tested  nu¬ 
merical  models.'’'6  Geographical  scales  for  ARAC 
assessments  vary  from  regional  (up  to  100  km)  to 
global  (thousands  of  km),  depending  on  the  re¬ 
lease  conditions.  At  present,  ARAC  services  are 
available  to  four  DOE  sites,  to  the  Federal  Avia¬ 
tion  Administration  (FAA),  to  DOE  emergency  re¬ 
sponse  operations  (e.g.,  the  Nuclear  Emergency 
Search  Team),  to  the  Nuclear  Regulatory  Com¬ 
mission's  Incident  Response  Center,  to  two 
nuclear  power  plants  (Indian  Point  and  Rancho 
Seco),  to  two  state  offices  for  emergency  services 
(New  York  and  California),  and  to  three  Depart¬ 
ment  of  Defense  (DOD)  sites. 

The  ARAC  center  (Fig.  2)  is  the  focal  point  for 
data  acquisition,  data  processing,  communica¬ 
tions,  and  assessments  through  the  use  of  inter¬ 
connected  minicomputers.  CDC  7600  computers 
at  the  LLNL  computer  center  are  used  to  calculate 
regional  and  global  atmospheric  transport  and  dif¬ 
fusion  estimates.  In  an  atmospheric  release  emer¬ 
gency,  the  ARAC  staff  can  obtain  exclusive  use  of 
a  CDC  7600  computer  within  a  matter  of  minutes 
after  the  computer  center  is  notified. 

Meteorological  data  from  the  Air  Force 
Global  Weather  Center  (AFGWC)  at  Offutt  Air 
Force  Base,  Bellevue,  Nebraska,  are  available  to 
the  ARAC  center  through  a  high-speed  data  link. 
The  ARAC  center  can  receive,  analyze,  display. 


and  store  meteorological  data  from  worldwide 
data  sources.  Selected  observational  and  forecast 
data  are  also  received  on  a  routine  scheduled  ba¬ 
sis.  In  an  emergency,  ARAC  can  request  supple¬ 
mental  data,  either  from  the  AFGWC  computer 
system  or  from  the  master  data  base  at  Carswell 
Air  Force  Base,  Fort  Worth,  Texas.  The  AFGWC 
gives  ARAC  high  priority  under  emergency  condi¬ 
tions,  thereby  speeding  up  the  response.  Cur¬ 
rently,  supplemental  and  backup  weather  data  are 
received  from  the  National  Weather  Service 
(NWS)  through  normal  teletype  and  facsimile 
channels.  In  the  future,  consideration  will  be 
given  to  including  the  updated  NWS  Automated 
Forecast  Office  Service  systems  in  the  ARAC 
network. 


ARAC  Models 

The  ARAC  system  utMfces  a  suite  of  numeri¬ 
cal  models  appropriate  to  a  variety  of  atmospheric 
release  incidents.  For  the  purpose  of  this  study, 
concentration  contours  coupled  with  the  SAI 
evacuation  model  were  calculated  by  using  the 
MATHEW  and  ADPIC  codes.  These  two  models 
are  described  briefly  below: 

MATHEW2** 

This  flexible,  three-dimensional  model  pro¬ 
vides  nondivergent  wind  fields  for  use  in  the 
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Figure  1.  ARAC  network. 
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Figure  2.  ARAC  Central  facility. 


ADPIC  model  and  in  other  studies.  It  adjusts  in¬ 
terpolated  wind-field  data,  subject  to  the  con¬ 
straints  of  mass  continuity  and  explicitly  intro¬ 
duced  topography.  Also,  it  is  available  for  other 
studies  such  as  assessing  windpower  sites  and  to¬ 
pographic  influences. 

rtUPIC4,5 

This  three-dimensional  particle-diffusion 
model  calculates  the  transport  and  diffusion  of  a 
puff  or  plume  in  a  time-varying  atmospheric 
boundary  layer.  It  is  based  on  the  particle-in-cell 
(PIC)  concept,  without  the  hydrodynamical  as¬ 


pects  of  the  conventional  PIC.  This  computer 
model  has  been  used  to  simulate  particulate  and 
gaseous  concentrations  from  one  or  more  sources 
at  distances  beyond  100  km,  general  deposition  of 
particles  with  given  size  distributions,  and 
rainout.  In  addition,  ADPIC  calculations  have 
been  compared  against  measurements  for  four 
different  field-diffusion  experiments.  The 
MATHEW/ADPIC  transport  and  diffusion  models 
are  continually  being  modified  and  verified 
against  field  tracer  studies  to  provide  ARAC  users 
with  useful  products  in  a  timely  manner. 


ARAC  Concentration  Contours 


The  ARAC  data  for  this  project  has  been  de¬ 
veloped  and  calculated  for  radionuclide  releases 
from  the  Rancho  Seco  Nuclear  Generating  Sta¬ 
tion.  Rancho  Seco  is  a  913-MWe,  pressurized  wa¬ 
ter  reactor,  nuclear  generating  station  operated 
since  1975  by  the  Sacramento  Municipal  Utilities 
District  (SMUD)  at  a  location  approximately  25  mi 
southeast  of  Sacramento,  California  (Fig.  3).  The 
site  is  somewhat  unusual  in  that  it  is  not  located 


adjacent  to  any  river  or  large  body  of  water,  and 
waste  heat  is  dissipated  to  the  atmosphere 
through  two  natural-draft  cooling  towers.  The 
area  surrounding  the  plant  site  is  almost  exclu¬ 
sively  agricultural.  More  densely  populated  areas 
of  greater  Sacramento  begin  about  10  mi  from  the 
site. 

The  plant  is  equipped  with  a  number  of  fixed 
radiation  monitoring  systems;  a  gamma-sensitive 


Figure  3.  Map  of  the  region  surrounding  the  Rancho  Seco  Nuclear  Power  Plant  C  —  25  mi 
southeast  of  Sacramento,  California). 


scintillation  detector  to  measure  radiation  levels  in 
the  spent  fuel  cooling  system;  high-energy 
gamma  detectors  to  monitor  gross  fuel  failure;  a 
gamma-sensitive  scintillation  detector  to  detect  ra¬ 


diation  levels  in  the  regenerant  holdup  tanks  dis¬ 
charge;  and  a  system  of  gross  beta  scintillation 
counters  and  single-channel  gamma  analyzers 
that  monitor  a  number  of  plant  processes,  the 


ventilation  systems,  and  the  atmosphere  both  in¬ 
side  and  outside  the  plant  tor  gaseous  and  air 
borne  particulate  radioactivity 

The  200 -ft  meteorological  tower  operated  b\ 
Rancho  Seco  is  located  about  3000  tt  east  of  tne 
reactor  building.  Wind  speed  and  direction  are 
measured  at  33  and  200  ft  above  ground  level 
Temperature  measurements  are  made  at  6  33  and 
200  ft  above  ground  level.  Wind  speed  and  direr 
tion  at  the  33- ft  level  are  measured  once  each  if' s 
the  other  parameters  are  measured  once  ea.  h 
minute. 

The  models  emploved  in  these  AKAC  cal¬ 
culations  require  the  topographv  on  the  region 
surrounding  the  Rancho  Seco  Nuclear  Power 
Plant  (NPP)  to  be  digitized  to  establish  a  site  map 
Figure  4  illustrates  the  topograph \  for  this  region 
in  three  dimensions.  The  isopleths  used  in  this  w 
port  were  computed  bv  running  MAI'HLU 


ADI ’ltd  li-r  .i  bi‘  km  grid  tor  a  release  time  of 
2  v-  ,■ 1  three  tv.a  i;»ii  s  Wile  ’chased.  '  Cs,  We. 
an  •  I  e  i,!.  ii’ii-jh'  i.ii-  ol  10  Ct/s.  10  Ci/s, 
and  i  t  i  •  r,  »i'eitivel’.  these  release  rates  were 
iho'ii.  bv  s\H.l)  as  r  presentative  ot  a  poten- 
:ial  i  cider.’  h.  till  event  of  ar  actual  nuclear  acci¬ 
dent  the  release  -ates  me. .  be  more  or  less.  The 
dose-  gner.  bv  ".e  AK  V.  calculated  contours 
would  then  attest  tlu  di-  .'ion  of  state  personnel 
to  ev .senate  oop- elation  ■.■■■tig  road  networks  that 
wt  uhi  mt-r-cct  cot. touts  vielding  dose  rates 
deemed  :oo  high  to  Oc-  su-.tained. 

I  igutes  A  ;  to  A  48  Appendix  A  give  the 
full  set  ot  contours.  The  externa!,  adult,  whole- 
bod  \  dose  spletb-  t  Xo  are  shown  in  Figs. 
A  1  to  A  It'  The  uhoh  boov  inhalation  dose 
contours  h-  C  -  are  ;  t  Ftgs  X  17  to  A-32.  The 
adult,  thvroia  rhalatioti  do--e  rate  contours  for 
I  are  in  1  :gs  \  33  to  A  4s 


Figure  4.  Three-dimensional  topography  of  the  region  surrounding  the  Rancho  Seco  Nuclear 
Power  Plant  (southwest  view).  Grid  origin  UTM  coordinates  are  x  -  646.0  km,  y  =  4237.0  km,  Z 
=  0  NASL.  Mesh  intervals  are  DELX  =  0.750  km,  DELY  =  0.750  km,  DELZ  -  25  m. 
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EVACD  Methodology 


Methodology  Overview 

A  rapid  computer  code  easily  adapted  for  use 
with  color  graphics  was  developed  to  convey  the 
dynamics  of  population  evacuation  in  the  vicinity 
of  a  nuclear  power  generating  station.  The  model 
developed  also  required  the  evacuation  model  to 
integrate  the  calculated  meteorological  transport 
and  diffusion  of  an  atmospheric  release  and  the 
resulting  dose  commitments  as  a  function  of  time 
and  position  from  the  release  point.  The  following 
section  discusses  the  methodology  developed  dur¬ 
ing  the  project. 

The  model  was  set  up  using  a  road  network 
consisting  of  links  and  nodes  overlaved  on  a  grid 
system  of  1600  cells  (40  X  40),  for  which  dose 
levels  were  calculated  for  each  time  step.  Geomet¬ 
rical  considerations  assured  the  road  network  and 
cells  were  coupled  together  correctly  for  all  cal¬ 
culations  requiring  both  reference  frames. 

Road  Network 

The  accommodation  of  vehicular  traffic  has 
been  the  primary  consideration  in  the  planning, 
design,  and  operation  of  streets  and  highways  by 
appropriate  groups  and  agencies.  Much  has  been 
written  and  compiled  about  the  requirements  nec¬ 
essary  to  successfully  design  and  implement  road 
networks  for  use  by  the  public.  This  subject  is 
broadly  referred  to  as  "highway  capacity,"  i.e.,  a 
measure  of  the  effectiveness  of  various  roadways 
to  accommodate  traffic.  The  determination  of 
highway  capacities  requires  a  general  knowledge 
of  traffic  behavior  and  a  specific  knowledge  of 
traffic  volumes  that  can  be  accommodated  under 
a  variety  of  roadway  configurations  and  operating 
conditions.  A  rational  and  practical  method  for 
determining  highway  capacities  has  always  been 
essential  for  an  economical  and  sound  utilization 
of  highway  transportation  systems.  Most  work 
performed  in  the  past  has  dealt  with  traffic  flows 
in  a  steady-state  environment.  To  model  a  dy¬ 
namic  evacuation  situation,  knowledge  gained 
from  the  steady  state  is  used  to  develop  a  time- 
dependent  traffic  flow  model. 

EVACD  Evacuation  Model 

General  Model  Overview 

EVACD  models  road  networks  as  a  collection 
of  roadways  and  intersections  that  can  be  repre¬ 


sented  by  links  and  nodes.  Figure  5  illustrates  a 
road  network  consisting  of  eight  nodes  and  nine 
links,  and  uses  the  nomenclature  of  the  EVACD 
model.  All  links  are  defined  by  the  nodes  that  are 
at  their  origin  and  termination.  The  subscripts  of 
the  links  identifier  /  are  first  ordered  as  origin 
node  number,  followed  by  the  termination  node. 
For  example,  a  link  extending  from  node  5  to  node 
4  is  identified  as  U4.  Travel  on  any  link  is  allowed 
in  only  one  direction.  Therefore,  roadways  with 
traffic  flow  in  both  directions  require  two  links,  as 
shown  by  and  in  Fig.  5.  The  arrows  in  Fig.  5 
indicate  the  direction  of  travel  on  each  link.  Thus, 
most  road  networks  can  be  easily  represented  as 
an  integrated  set  of  nodes  and  links. 

A  unique  feature  of  the  link  network  is  that 
each  is  defined  by  a  set  of  descriptors  that  gives 
all  data  parameters  particular  to  each  link.  The 
descriptors  include  the  following: 

1.  Origin  node. 

2.  Termination  node. 

3.  Number  of  travel  lanes. 

4.  Length  of  link. 

3.  Free  flow  speed. 

6.  Standard  traffic  capacity. 

7.  Jam  density. 

8.  All  factors  affecting  traffic  flow. 

The  model  in  EVACD  does  not  track  individ¬ 
ual  vehicles,  but  rather  uses  relationships  between 
flows,  road  characteristics  such  as  speed,  density, 
and  traffic  backup,  and  other  relevant  traffic  pa¬ 
rameters.  EVACD  is  not  intended  to  model 
steady-state,  random  traffic  problems.  It  was  de¬ 
veloped  explicitly  to  solve  complex  evacuation 
problems.  The  model  is  based  on  the  concepts 
found  in  Retf.  1. 


Figure  5.  Representation  of  a  road  network 
using  the  nomenclature  of  the  EVACD  model. 
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FVACD  has  been  designed  to  display  the  dy¬ 
namics  of  each  portion  of  a  traffic  network  at 
snapshots  in  time  during  an  evacuation.  The  dy¬ 
namics  found  on  each  link  and  at  each  node  are 
easilv  displaved.  enabling  identification  of  possi¬ 
ble  problem  areas  during  an  actual  evacuation. 
EVACD  has  also  been  designed  so  that  modifica¬ 
tions  and  additions  can  easily  be  made  in  the  fu¬ 
ture.  (See  Appendix  B  for  a  detailed  discussion  of 
the  EVACD  computer  program  and  Appendix  C 
for  a  listing  of  the  EVACD  main  routine,  subrou¬ 
tines  and  common  tiles  ) 

EVACD  Model  Dynamics 

EVACD  operates  as  a  numerical  approxima¬ 
tion  scheme  to  the  total  evacuation  process.  By 
combining  the  relationships  of  the  road  network 
parameters  to  a  mass  balance  of  vehicles  on  the 
network  over  short  periods  of  time,  an  approxi¬ 
mation  of  the  evacuation  can  be  calculated  in  a 
step-wise  manner.  This  is  accomplished  first  by 
determining  the  dynamics  on  each  link  during  a 
short  period  of  time  or  during  a  time  step.  The 
conditions  on  the  link  at  the  end  of  each  time  step 
are  calculated,  such  as  new  road  densities  and 
flow  speeds  and  the  number  of  vehicles  traversing 
and  reaching  the  end  of  the  links.  Afterwards,  a 
mass  or  vehicular  balance  is  performed  at  the 
nodes,  and  the  number  of  vehicles  moving 
through  each  node  from  one  link  to  another  is 
calculated.  Thus,  a  time  step  typically  involves  the 
calculation  of  the  dynamics  on  all  links  followed 
by  a  vehicular  balance  at  each  node.  If  the  time 
steps  are  chosen  appropriately,  the  evacuation 
process  can  be  modelled  over  time  by  stepping 
through  each  time  step  individually. 

As  can  be  seen,  the  model  essentially  in¬ 
volves  two  calculational  procedures:  a  dynamical 
calculation  for  all  links  during  a  time  step  and  a 
vehicular  balance  at  all  nodes  at  the  end  of  a  time 
step.  These  procedures  are  identified  as  the  dy¬ 
namic  link  scan  and  the  vehicular  balance  node 
scan. 

Dynamic  Link  Scan.  During  the  link  scan,  the 
traffic  volume  is  calculated  from  the  beginning 
condition  on  each  link.  The  network  shown  in 
Fig.  5  illustrates  each  link.  A  link  can  have  one  or 
more  lanes  (defined  by  NL,  )  where  lmn  is  the 

link  designator  as  discussed  for  Fig.  5.  To  simplify 
writing  the  appropriate  relations  and  equations, 
the  nomenclature  will  be  link  /  instead  of  /mn,  in 
which  l  represents  a  number  corresponding  from 
one  to  the  total  number  of  links,  n,. 


Thus,  the  number  of  traffic  lanes  for  link  /  is 
NL,.  The  length  of  link  I  in  miles  is  LD:.  The  free¬ 
flow  operating  speed  on  link  I  is  FFS,,  given  in 
miles  per  hour.  The  free-flow  operating  speed  is 
the  operating  speed  of  vehicles  on  the  link  during 
extremely  low  traffic  densities.  The  standard  ca¬ 
pacity  of  link  /  is  CS(/),  given  in  vehicles  per  hour. 
Standard  capacity  is  the  maximum  number  of  ve¬ 
hicles  that  can  reasonably  be  expected  to  pass 
over  a  given  section  of  roadway  in  1  h.  (These  and 
other  parameters  are  discussed  thoroughly  in 
Ref.  2.) 

To  expand  on  the  model,  a  few  general  traffic 
flow  relationships  must  first  be  explained.  The 
three  key  link  parameters  are 
F  =  rate  of  flow  (vehicles/h), 

S  =  average  link  speed  (mi/h),  and 
D  =  link  density  per  lane  (vehicles/mi/lane). 

The  parameters  are  related  via  the  expression 

F  =  S  ■  D  , 

which  essentially  says  the  number  of  vehicles 
traveling  on  a  link  is  equal  to  the  average  speed 
multiplied  by  the  link  density.  These  parameters 
are  macroscopic  measures  of  a  traffic  flow.  Al¬ 
though  the  previous  expression  seems  to  suggest 
that  a  given  rate  of  flow  may  occur  at  numerous 
combinations  of  speed  and  density,  this  is  not  the 
case.  In  practice,  only  a  limited  number  of  com¬ 
binations  will  occur  since  there  may  be  additional 
relationships  between  F  and  S,  F  and  D,  and  S  and 
D  which  control  these  combinations.  Figure  6 
graphically  illustrates  the  general  form  of  these 
three  relationships.  (References  7  and  8  give  an  in- 
depth  discussion  of  these  relationships.) 

During  the  link  scan,  the  density  on  each  link 
is  calculated  first,  where  t  is  equal  to  time.  For 
illustration,  t  is  taken  to  represent  time  at  the  end 
of  a  time  step.  The  link  density  per  lane  is 

L,  ( t)/NL, 

(t>  LD,  -  LQ(t)  ' 
where 

V,(t)  =  number  of  vehicles  moving  on  link  /, 
and 

LQ,(t)  =  length  of  the  queue  of  vehicles  stopped 
at  the  end  of  link  l  because  of  a  bottle¬ 
neck  (mi). 

IQ,  is  the  result  of  traffic  bottlenecks  due  to 
flow  restrictions  at  the  end  of  a  link.  For  example, 
two  links  may  end  at  a  node  with  only  one  exit 


Figure  6.  Relationship  between  speed,  flow, 
and  density.78  Flow  rate  F,  occurs  under  two 
different  flow  conditions,  A  and  B. 


link.  Thus,  vehicles  may  be  required  to  stop  at  the 
ends  of  the  input  links  because  of  a  flow  restric¬ 
tion.  LQ i  is  the  length  of  such  a  queue  of  vehicles 
at  the  end  of  link  /.  If  VQ,(t)  is  the  number  of 
vehicles  in  the  queue  and  VL  is  vehicle  length  in 
miles,  then 

LQ,(t)  =  VQ,(t)  •  VL  . 

The  average  link  speed  is  calculated  in  the 
EVACD  model  by  assuming  a  parabolic  relation¬ 
ship  between  the  link  flow  and  density,  shown  in 
Fig.  6.  The  jam  density  per  lane  D],  illustrated  in 
Fig.  6  can  be  approximated  as 


Therefore,  the  average  link  speed  is  defined  as 


S,(t)  =  FFS,  1 1  - 


and  the  link  rate  of  flow  becomes 


F,(t)  =  D,(t)  ■  S,(t)  ■  NL,  . 

The  number  of  vehicles  reaching  the  termination 
node  or  the  end  of  the  queue  during  the  interval  is 
defined  as 


VL,(t)  =  F{t)  ■  T  , 

where  T  =  time  step  interval  (h). 

VL^t)  is  constrained  to  be  equal  to  or  less  than  the 
number  of  vehicles  initially  on  the  link,  or 

VL,(t)  <  V,(t)  . 

The  vehicles  remaining  on  the  link,  VR,,  are  de¬ 
fined  to  be 

VR,(t)  =  V,(t)  -  VL,(t)  . 

After  the  time  step,  the  link  may  still  have  room 
to  add  additional  vehicles  during  the  next  scan. 
This  excess  vehicular  capacity  VE,  is  calculated  by 
assuming  that  the  added  vehicles  will  essentially 
bring  the  link  density  D;,  to  the  jam  density  D/,, 

VE,  (f  i 

{W^LQ^\Tnl,  ' 

or 

VE, (ft  =  [LD,  -  LQ, ifi]  •  NL,\DJ,  -  Dtt  \  . 

The  vehicular  balance  node  scan  follows  comple¬ 
tion  of  the  dynamic  link  scan. 

Vehicular  Balance  Node  Scan.  The  vehicular 
node  scan  calculates  the  number  of  vehicles  flow¬ 
ing  in  and  out  of  each  node,  assuring  a  vehicular 
balance  over  time.  The  process  starts  by  summing 
up  all  potential  vehicles  that  may  use  a  node, 

VNjfl  =  VQ/fi  +  VI/ f  l  , 

where  n,  =  subscript  for  node  n  flow  from  link  /. 

Thus,  VN,,^t)  is  the  sum  of  the  end  queue  of  link  / 
plus  all  vehicles  reaching  the  input  link  queue  or 
the  link  termination  during  the  time  step. 

Next,  it  is  necessary  to  calculate  the  fraction 
of  time  that  traffic  flow  from  an  input  link  to  the 
node  can  move  through  the  node  or  intersection, 
G,.  For  cases  such  as  signalized  intersections,  the 
value  may  be  defined  as 

C,(0  =  GS,  , 

where  GS,  =  fraction  of  time  for  flow  (fraction 
h/h). 
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For  cases  where  there  is  no  signalized  intersection, 
the  fraction  of  flow  time  for  a  link  is  assumed  to 
be  proportional  to  its  potential  flow  to  those  of  all 
other  input  links  to  node  n, 

„  VNJD/NL, 

wjl  *  * - 

V  VNnk/NLk 


where  k  =  all  input  links  to  node  n. 

The  approach  capacity  of  each  input  link,  CA,,  is 
calculated  to  be 


Next,  it  is  necessary  to  determine  the  number 
of  vehicles  passing  through  the  node  into  each 
output  link.  The  calculational  procedure  is  as  fol¬ 
lows.  First,  the  number  of  vehicles  leaving  each 
input  link  is  calculated  assuming  no  node 
restraints. 

Vim  «  T  ■  CA,(t) 

and  is  restricted  to  be  equal  to  or  less  than  VN 
the  number  of  vehicles  available  for  leaving  the 
link.  Thus, 

VIM  =  VNJt),  if  T  •  CA,(t)  <  VNJt)  . 

The  flow  volume  received  by  each  output 
link  from  node  n  is  ascertained  by  first  calculating 
a  preference  for  traffic  to  take  the  link.  We  assume 
that  the  preference  is  proportional  to  a  previously 
agreed  upon  preference  factor  multiplied  by  the 
average  link  speed. 


PF,  •  S/f! 


p,m  -  ■" 

^PF*  •  Sk\l) 


where 

P,m  —  preference  factor  from  input  link  /  to 
output  link  j  (fraction), 

PF,  =  inputted  preference  factor  for  output 
link  /  regardless  of  input  link,  and 
k  —  all  output  links  to  node  n. 

The  vehicles  received  by  output  link  VO,  are 
defined  as 


VO/ fl  =  y  Vlt(/l  •  PA/tl  , 


where  k  =  all  input  links  to  node  n. 

This  number  must  be  less  than  or  equal  to  the 
capacity  of  output  link  j  multiplied  by  time  step  T 
and  cannot  exceed  VE,(f),  the  excess  vehicular  ca¬ 
pacity  on  the  link.  This  can  be  expressed  as 

lO,(t)  =  min  [VO,(ty,  T  ■  CS VE,<f)]  , 


where  min 
quantities. 


take  the  minimum  of  all  the 


Therefore,  the  vehicles  VTU  transferred  from  input 
link  /  to  output  link  j  are 


vT,m  =  vi,(t> .  p(/t  i 


y  vihw  ■  m 


where  k  =»  all  input  links  to  node  n. 

The  vehicular  node  scan  is  complete  follow¬ 
ing  this  calculation.  The  process  continues  with 
another  dynamic  link  scan  for  the  next  time  step 
where 


V/fl  =  VRfit  -  Ti  +  y VTkft  -  Ti  , 

4 7 * 

where  k  =  all  input  links  to  origin  node  of  link  /. 

Therefore,  by  time-stepping  through  the  evacua¬ 
tion  and  performing  a  dynamic  link  scan  followed 
by  a  vehicular  balance  node  scan  for  each  time 
step,  the  desired  time-dependent  results  can  be 
calculated. 


Calculation  of  Appropriate  Doses 

The  EVACD  model  described  in  this  report 
calculates  three  specific  radiological  components 
for  assessing  the  effect  of  airborne  radioactive  re¬ 
leases  from  a  nuclear  power  generating  station. 
The  calculations  are  contained  in  subroutine 
DSCLC.  The  calculations  were  performed  for  a 
specific  set  of  data: 

1.  Accumulated  dose  commitments  (mrem) 
to  the  thyroid  due  to  13II  inhalation. 
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2.  Accumulated  dose  commitments  (mrem) 
to  the  total  body  due  to  U7Cs  inhalation. 

3.  Accumulated  whole-body  dose  exposures 
(mrem)  to  13,Xe. 

For  other  applications,  DSCALC  could  be 
modified. 

The  calculated  population  doses  to  the  evac¬ 
uating  population  were  totaled  for  each  time  step 
by  adding  up  the  following  components: 


1.  All  populations  at  input  or  output  nodes. 

2.  All  populations  at  link  output  queues 

3.  All  populations  traveling  on  links  during 
the  time  step. 

Geometric  considerations  in  the  subroutines 
LNKSET  and  NODSET  assure  calculation  of  the 
correct  dose  additions  for  each  link  and  each 
node. 


Sample  Case  and  Outputs 


The  evacuation  emergency  response  model 
coupling  ARAC  output  was  run  for  a  test  case  for 
the  region  surrounding  the  Rancho  Seco  NPP,  lo¬ 
cated  ~~25  mi  southeast  of  Sacramento.  The  case 
was  set  up  with  three  bottlenecks  to  restrict  the 
flow  capacity  of  certain  links.  The  population  of 
the  region  within  a  10-mi  radius  of  the  nuclear 
power  plant  is  not  dense  enough  to  test  the  ca¬ 
pabilities  of  this  model.  That  is,  for  the  given 
population  and  roads  surrounding  Rancho  Seco, 
the  population  would  be  evacuated  within  about 
1  h.  A  sample  problem  employing  evacuation  of  a 
population  of  10,304  people  on  roads  with  limited 
traffic  capacities  is  described  in  Appendices  D  and 
E,  which  detail  model  parameters,  population  dis¬ 
tribution,  and  output.  These  data  have  enabled  us 
to  fully  display  this  evacuation  emergency 
sponse  model  on  computer  graphics  terminals. 

The  EVACD  program  was  developed  for  use 
on  a  VAX-780  in  conjunction  with  a  VSV11  color 
graphics  terminal  made  by  Digital  Equipment 
Corp.  The  graphics  package  used  was 
GRAFCORE,9  developed  by  the  Lawrence  Liver¬ 
more  National  Laboratory.  The  output  of  the  pro¬ 
gram  consists  of  three  types  of  displays.  The  first 
type  appears  only  once  and  shows  the  evacuation 
zone  as  a  road  network.  The  second  type  shows 
the  road  network  with  the  isopleths  the  user 
wishes  to  see.  The  third  type  is  a  histogram  show¬ 
ing  population  movement. 

In  the  displays  thick  lines  denoting  queues 
build  up  behind  the  nodes.  During  later  time  steps 
these  same  queues  gradually  shrink.  The  displays 
show  the  growth  and  dispersion  of  the  radioactive 
plume  (represented  by  the  dose  isopleths).  The 
user  can  see  which  nodes  and  links  intersect  the 
path  of  the  plume  and  how  the  population  can  be 
diverted.  Growth  and  shrinkage  of  nodes  and 
links  provide  information  on  the  movement  of 
people  and  on  the  size  of  the  population  at  risk 
from  the  plume.  The  histograms  show  population 
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movement  outward  and  its  slow  movement 
through  the  rings  where  the  bottlenecks  are. 

Figure  7  is  a  black-and-white  illustration  of 
the  first  type  of  display  showing  the  initial  road 
network,  site  layout,  and  evacuation  case  param¬ 
eters  before  evacuation  has  proceeded  It  is  the 
first  display  available  to  the  user  of  this  model. 
The  reactor  building  and  site  boundary  are  dis¬ 
played  in  orange  in  the  center  of  the  screen;  a  lake 
within  the  site  is  cyan.  The  axes  show  the  distance 
in  mi  from  the  reactor  building.  Two  circles  in 
green  are  at  radii  of  10  and  13  mi,  respectively. 
The  road  network  is  displayed  in  yellow  and  the 
small  circles  or  nodes  are  in  orange.  The  road  net¬ 
work  ends  beyond  the  10-mi  radius,  at  which  dis 
tance  the  population  is  assumed  to  be  evacuated. 

At  the  beginning  of  the  problem,  all  road 
thicknesses  and  node  circles  are  the  same  size.  As 
the  evacuation  proceeds,  the  thickness  of  the  road 
may  increase  indicating  a  queuing,  and  the  circles 
increase  in  size  indicating  a  bottleneck  of  the 
population  at  these  nodes.  In  the  upper  left-hand 
portion  of  the  screen  are  the  warning  time  of 
0.25  h  required  to  evacuate,  the  population  time  of 
0  h  needed  to  evacuate  after  the  warning  is  given, 
and  the  total  population  contained  within  a  10-mi 
radius  of  the  plant.  The  user,  of  course,  has  the 
option  of  varying  the  warning  time  and  the  prepa¬ 
ration  time  if  other  numbers  are  more  appropriate 
to  an  exercise  or  real  accident. 

In  the  second  type  of  display,  the  user  is 
prompted  for  the  particular  nuclide  whose  iso¬ 
pleths  he  wishes  to  see.  These  isopleths  are  dis¬ 
played  over  the  existing  road  netw  ork.  The  user  is 
then  asked  if  he  wants  to  see  isopleths  for  a  differ¬ 
ent  nuclide.  If  the  answer  is  yes,  the  user  chooses 
the  next  nuclide  to  be  displayed  and  the  previous 
plot  is  replaced.  If  the  answer  is  no,  the  display 
disappears  and  the  histogram  plot  is  displayed. 

The  second  type  of  display  is  similar  to  the 
first  but  changes  dynamically  at  each  dispersion 
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Figure  7.  First  type  of  display  of  the  EVACD  program  show* 
ing  initial  road  network,  site  layout,  and  evacuation  case 
parameters. 


update  time  as  the  evacuation  proceeds.  Orange 
circles  are  drawn  around  each  input  node.  Gold 
circles  are  drawn  around  each  output  node  to  in¬ 
dicate  the  number  of  people  evacuated.  Links  are 
redrawn  as  thicker  lines  to  indicate  road  density, 
and  queues  are  drawn  in  orange.  In  addition,  the 
user  can  request  that  dose  isopleths  be  added  to 
the  display. 

Figure  8,  the  second  type  of  display  (road 
network-isopleth),  shows  dose  isopleths  due  to 
,31I  inhalation  0.25  h  following  release.  The  pro¬ 
gram  counts  zero  people  on  the  road  network 
since  the  population  has  not  yet  begun  to  move. 
Four  isopleths  are  displayed,  corresponding  to  the 
four  contour  levels  read  from  the  input  file 
ACLEVS.DAT.  These  levels  are  typed  in  the  top 
right  comer  of  the  screen.  On  the  color  graphics 
terminal  the  isopleths  are  drawn  in  varying 
shades  of  magenta,  with  the  darkest  shade  indi¬ 
cating  the  highest  dose.  The  values  of  the  levels 
shown  on  the  screen  are  also  typed  in  shades  of 
magenta  corresponding  to  the  appropriate 
isopleth. 


Further  information  in  the  top  left  comer  of 
the  display  includes  time  since  release,  the  popu¬ 
lation  entering  the  road  network  at  the  time  of 
evacuation,  the  population  leaving  the  road  net¬ 
work  at  that  time,  and  the  total  fraction  of  the 
population  evacuated  since  the  release. 

The  third  type  of  display  occurring  at  each 
dispersion  update  is  a  population  histogram 
showing  the  number  of  people  presently  in  10 
rings  around  the  site.  The  rings  are  1-mi  wide,  i.e., 
from  0-1  mi,  1-2  mi,  9-10  mi.  A  bar  showing  the 
number  of  people  evacuated  beyond  10  mi  is  also 
displayed.  Figure  9  (a  and  b)  is  an  example  of  the 
second  and  third  types  of  displays,  where  Fig.  9a 
shows  the  desired  isopleths  at  0.50  h,  and  Fig.  9b 
is  the  corresponding  histogram.  In  addition  to  the 
time  and  population  information  in  the  top  left 
comer  of  the  histogram,  accumulated  population 
doses  for  131I,  133Xe,  and  137Cs  appear  in  the  top 
right  comer.  These  values  are  calculated  as  de¬ 
scribed  in  "Subroutine  DSCALC." 

Figures  10  through  21  represent  a  black-and- 
white  version  of  the  color  graphics  output  from  a 
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Figure  8.  Second  type  of  display  of  the  EVACD  program 
showing  l3lI  inhalation  dose  isopleths  at  0.25  h  (15  min)  fol¬ 
lowing  release.  In  this  type  of  display  the  user  requests  spe¬ 
cific  isopleths  to  be  shown. 


sample  execution  of  the  EVACD  program,  show¬ 
ing  in  sequence  (in  sets  of  road  network-isopleth 
and  corresponding  histogram  displays)  the  evac¬ 
uation  of  the  population  through  the  road  net¬ 
work  every  15  min  up  until  3  h  after  the  release, 
when  85%  of  the  population  is  evacuated.  Iso- 
pleth  displays  for  13  I  were  requested  at  every  dis¬ 
persion  update,  while  isopleth  displays  for  l33Xe 
and  l37Cs  were  requested  at  every  third  dispersion 
update.  The  user  of  the  program  has  the  option  of 
skipping  any  of  these  sets  or  backtracking  to  help 
in  the  decision  making.  At  each  display  of  the 
road  network,  the  user  has  the  option  of  deploy¬ 
ing  the  radiation  isopleths  for  131 1,  133Xe,  137Cs,  or 
none  at  all. 

The  first  set  of  displays  (road  network- 
isopleth  and  corresponding  histogram)  in  Fig.  10 
shows  the  evacuation  0.75  h  after  the  evacuation 
has  started.  In  the  central  portion  of  the  road  net¬ 
work-isopleth  display,  isopleths  of  137Cs  inhala¬ 
tion  dose  rates  (mrem)  extend  from  the  inner  con¬ 
tour  of  0.01  mrem  to  the  outer  one  of  1  X  10  3 


mrem.  These  contours  spread  to  the  upper  left- 
hand  portion  of  the  display,  intersecting  person¬ 
nel  evacuated  along  these  routes.  Thus,  emer¬ 
gency  response  personnel  may  decide  to  change 
their  evacuation  plans  based  on  these  factors.  The 
circles  outside  the  10-mi  radius  have  expanded 
proportionally  to  the  population  numbers  evacu¬ 
ated  to  that  node.  Data  in  the  upper  left-hand  cor¬ 
ner  inform  the  user  that  1042  people  are  still  en¬ 
tering  the  road  network,  6808  people  have  left  the 
road  network,  and  a  total  of  66.1%  have  already 
been  evacuated  outside  the  10-mi  radius. 

Note  the  large  circles  surrounding  some  of 
the  nodes  in  the  left-central  portion  of  the  display. 
A  thickening  of  the  road  network  is  also  visible 
along  the  lower  left-central  and  upper  left-central 
portion  of  the  graph.  The  large  circles  and  the 
thickening  relate  to  the  corresponding  histogram 
which  shows  the  number  of  people  still  being 
evacuated  as  a  function  of  distance  from  the  plant 
and  those  already  evacuated  45  min  after  evacua¬ 
tion  began.  The  bottlenecks  between  8  and  10  mi 
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Figure  9.  Set  of  (a)  road  network-isopleth  and  (b)  corresponding  histogram  displays  at  0.50  h.  The 
histogram,  depicting  population  movement,  is  the  third  type  of  display  of  the  EVACD  program. 
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Figure  10.  (a)  Road  network* 
isopleth  display  and  (b)  cor¬ 
responding  histogram  of  a  test 
case  for  the  region  surrounding 
the  Rancho  Seco  Nuclear  Power 
Plant  showing  evacuation  of  the 
population  through  the  road  net¬ 
work  at  0.75  h.  Figures  11 
through  21  continue  the  se¬ 
quence  showing  evacuation  ev¬ 
ery  15  min  up  until  3  h  after  the 
release,  when  85%  of  the  popula¬ 
tion  is  evacuated.  This  sequence 
is  a  black-and-white  version  of 
the  color  graphics  output  from  a 
sample  execution  of  the  EVACD 
program. 
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and  3  and  5  mi  can  be  seen  now  on  the  histogram 
and  are  related  to  the  circles  seen  in  the  cor¬ 
responding  road  network-isopleth  display.  The 
number  evacuated  is  shown  in  the  upper  left- 
hand  comer  of  the  histogram  and  is  directly  re¬ 
lated  to  the  circle  outside  the  10-mi  radius.  Total 
accumulated  population  doses  are  displayed  in 
the  upper  right-hand  comer. 

In  this  particular  case,  the  Rancho  Seco  test 
case,  the  sequence  is  continued  until  the  program 
terminates  with  92.6%  of  the  population  evacu¬ 
ated  (see  the  final  set  of  displays  in  Fig.  21).  The 
user  has  the  option  of  stopping  this  program  at  a 
preset  particular  percentage  of  population  evacu¬ 
ated  or  at  any  time. 

In  summary,  information  is  available  to  the 
user  as  a  function  of  time  to  assist  in  either  train¬ 
ing  sessions  or  an  actual  emergency  response.  At 
the  beginning  of  an  accident  and  for  every  15  min 
of  simulated  time  thereafter,  the  user  may  do  the 
following: 

•  View  a  colored  map  indicating  the  road 
network  and  the  number  of  individuals 


within  a  10-  and  15-mi  radius  of  the  acci¬ 
dent  at  the  NPP. 

Display  dose  contours  for  U1l,  B3Xe,  and 
B/Cs  overlayed  on  the  colored  map. 
Display  a  histogram  showing  the  number 
of  people  still  within  the  evacuation  zone, 
those  evacuated,  and  accumulated  popu¬ 
lation  doses. 

Determine  directly  from  the  color  graphics 
displayed  whether  a  bottleneck  exists  in 
evacuating  individuals  along  a  particular 
clogged  road. 

Determine  directly  from  the  color  graphics 
and  overlayed  contours  whether  a  par¬ 
ticular  evacuation  route(s)  intersects  the 
radiation  plume  endangering  individuals 
being  evacuated. 

Skip  any  of  the  preceding  options  or 
backtrack  to  further  study  a  particular 
display. 


Suggested  Improvements 


The  EVACD  model  was  developed  to  dem¬ 
onstrate  the  feasibility  of  coupling  an  evacuation 
model  with  a  color  graphics  system  to  illustrate 
the  useful  information  available  from  such  a  com¬ 
bination.  The  display  of  isopleths  resulting  from 
an  airborne  radioactive  release  using  ARAC  data 
to  show  potential  evacuation  problems  was  an  in¬ 
tegral  part  of  the  demonstration.  Coupling  of  the 
evacuation  model  with  the  airborne  releases  was 
demonstrated  by  calculating  population  doses  and 
dose  commitments.  The  objective  of  the  program 
was  achieved  and  has  enabled  improvement  of 
the  system.  The  following  are  suggestions  for  im¬ 
provement  to  the  developed  model: 

1.  Add  more  interactive  capabilities  to  allow 
operators  running  the  model  to  modify 
road  network  conditions  at  each  appro¬ 
priate  time  step.  In  otheT  words,  the  oper¬ 


ator  should  be  allowed  to  insert  road 
blocks,  change  road  preference  factors, 
and  similar  types  of  real-time  changes. 
Add  more  graphics  to  make  the  model 
the  basic  tool  for  training  emergency  pro 
cedure  personnel.  Suggested  output 
would  aid  in  assessing  the  various  evac¬ 
uation  decisions. 

Generalize  the  dose  assessment  routines 
to  cover  more  cases  of  a  general  nature. 
Refine  the  road  network  model  to  include 
more  specific  data,  such  as  road- 
dependent  capacities  similar  to  those 
found  in  Ref.  10. 

Calibrate  the  model  against  data  for  past 
actual  evacuations. 

Develop  a  training  program  using 
EVACD  as  the  main  teaching  tool. 
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Appendix  A: 

ARAC  Contours 

Appendix  A  contains  the  full  set  of  ARAC  calculated  contours  for  the  external,  adult  whole-body 
dose  isopleths  for  133Xe  (Figs.  A-l  through  A-16);  the  whole-body,  inhalation  dose  for  13  Cs  (Figs.  A- 17 
through  A-32);  and  the  adult,  thyroid,  inhalation  dose  rate  for  !3,I  (Figs.  A-33  through  A-48). 
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Figure  A-21. 


Figure  A-23. 
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Figure  A-24, 


Figure  A-27 


Figure  A-26. 


Figure  A-28. 
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Figure  A-30. 


Figure  A-32. 
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Figure  A-42. 


Figure  A-44. 
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Figure  A-45. 
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Appendix  B: 

Discussion  of  EVACD  Computer  Program 

Appendix  B  provides  a  description  of  the  EVACD  computer  program  that  simulates  the  methodology 
discussed  earlier  in  "ARAC."  The  first  part  of  Appendix  B  describes  the  main  routine  called  EVACD, 
which  is  followed  by  an  alphabetical  listing  of  each  subroutine.  Some  of  the  subroutines  make  reference 
to  color  graphics  displays;  these  displays  are  reported  in  detail  in  the  text  in  "Sample  Case  and  Outputs." 
Each  common  block  of  the  program  appears  in  a  separate  file;  INCLUDE  statements  were  used  to  incor¬ 
porate  these  common  blocks  in  the  appropriate  routines.  The  second  part  of  Appendix  B  describes  the 
input  files  required  by  the  program.  A  test  case  using  the  sample  input  files  was  run;  the  output  is 
described  in  the  text  in  "EVACD  Methodology." 


Description  of  EVACD  Main  Routine  and  Subroutines 

Main  Routine  (EVACD) 

EVACD  is  the  main  driver  routine  for  the  program.  It  is  responsible  for  calling  each  subroutine  in  the 
correct  order.  In  addition,  much  of  the  calculations  and  reading  of  data  occurs  in  EVACD.  Specifically, 
data  describing  the  road  network  of  links  and  nodes,  plus  initial  population  placement,  are  read  from  the 
file  EVACDB.INP  on  unit  10.  Information  describing  the  specific  evacuation  case  to  be  run  is  read  from 
the  file  EVACDB.CAS  on  unit  11.  This  file  includes  data  which  describe  the  evacuation  zone — a  15-  by 
15-mi  square-shaped  zone  divided  into  1600  grid  squares  in  a  40  X  40  arrangement.  EVACDB.CAS 
contains  the  x-  and  v-coordinates  of  the  lower  left-hand  comer  of  the  grid  and  the  x-  and  y-dimensions  of 
each  grid  square. 

EVACD  is  arranged  in  four  main  sections.  The  first  section,  which  is  executed  only  once,  reads  the 
data  from  units  10  and  11,  calculates  the  location  of  each  node  on  the  40  X  40  grid  and  the  length  of  each 
link  and  calculates  the  total  population  to  be  evacuated.  Then  in  the  second  section,  EVACD  enters  a  time 
increment  loop,  which  is  executed  a  number  of  times  as  specified  in  the  EVACDB.INP  input.  Within  this 
time  loop  are  calls  to  the  various  subroutines  that  provide  the  dynamic  color  graphics  display  of  the 
evacuation.  The  dynamic  display  is  not  updated  at  every  time  step;  it  is  updated  only  once  within  a  time 
interval  called  "dispersion  update  time,"  specified  in  EVACEDB.CAS.  Currently  the  update  occurs  every 
15  min. 

The  third  and  fourth  sections  are  a  link  and  a  node  loop,  both  of  which  are  nested  within  the  time 
loop.  The  link  loop  performs  the  dynamic  link  scan  as  described  in  "EVACD  Methodology:  Dynamic  Link 
Scan."  Variables  calculated  include  the  population,  flow,  and  queue  associated  with  each  link  at  each  time 
step.  The  node  loop  performs  the  vehicular  balance  node  scan  described  in  "EVACD  Methodology: 
Vehicular  Balance  Node  Scan."  Variables  calculated  include  the  population  entering  and  leaving  each 
node  at  each  time  step. 

The  time  loop  is  terminated  when  all  of  the  numbers  of  specified  time  steps  have  been  executed,  or 
when  a  specified  fraction  of  the  population  has  been  evacuated  beyond  a  10-mi  radius  from  the  release 
point.  This  fraction  is  also  input  from  EVACDB.CAS.  The  program  ends  following  termination  of  the  time 
loop. 

Subroutine  DATA 

At  each  dispersion  update  (15  min  intervals)  DATA  is  called.  This  subroutine  reads  the  airborne 
release  meteorological  transport  and  diffusion  data  supplied  by  Lawrence  Livermore  National  Laboratory. 
The  data  are  presented  in  three  separate  input  files.  File  1131. LLNL  on  Unit  15  contains  thyroid  dose  due 
to  inhalation  of  l3,I.  File  XE133.LLNL  on  unit  16  contains  whole-body  doses  due  to  exposure  to  133Xe.  File 
CS137.LLNL  on  unit  17  contains  total-body  doses  due  to  inhalation  of  137Cs.  In  each  case,  the  doses  are 
given  in  mrem  at  each  of  the  1600  grid  squares.  The  doses  change  every  15  min  until  4.5  h  after  the  release. 

Subroutine  DISP 

DISP  is  called  at  the  start  of  the  calculation  and  at  each  dispersion  update.  It  displays  the  current  link 
and  node  numerical  data.  The  size  of  the  population  density  on  each  link  is  reflected  by  drawing  the  link 
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as  a  thicker  line  for  higher  densities.  Similarly,  dense  queues  are  drawn  as  thick  lines.  On  the  color 
graphics  display,  links  are  drawn  in  yellow  and  queues  are  drawn  in  orange.  The  population  entering  each 
input  and  output  node  is  symbolized  by  four  equally  spaced  concentric  circles.  The  radius  of  the  outer 
circle  reflects  the  size  of  the  population  of  the  node,  with  larger  radii  indicating  larger  populations.  Input 
nodes  are  drawn  in  orange  and  output  nodes  are  drawn  in  gold. 

Subroutine  DSCALC 

DSCALC  is  called  at  each  step  following  the  execution  of  the  link  loop.  DSCALC  calculates  the  dose 
from  1311,  133Xe,  and  137Cs  to  the  population  presently  on  each  link.  This  is  done  by  comparing  the  grid 
boxes,  through  which  a  particular  link  passes,  with  the  doses  at  the  center  of  each  box  (input  files 
1131.LLNL,  XE133.LLNL,  and  CS137.LLNL).  At  each  time  step  the  total  accumulated  doses  from  the  three 
nuclides  are  calculated.  At  each  dispersion  update  these  total  doses  are  typed  on  the  color  graphics  display 
in  the  top  right  comer. 

Subroutine  ERROR 

Some  of  the  information  required  by  the  program  is  provided  interactively  by  the  user.  If  the  user 
makes  a  mistake,  ERROR  will  type  an  error  message  and  allow  the  user  to  try  again. 

Subroutine  HISTO 

HISTO  provides  a  color  graphics  display  of  a  histogram  plot  depicting  the  population  in  1-mi  rings 
around  the  reactor  site  at  each  dispersion  update.  The  rings  go  out  to  10  mi,  with  an  eleventh  histogram 
showing  the  number  of  people  evacuated  beyond  the  10-mi  zone.  This  plot  is  useful  because  it  gives  an 
indication  of  the  progress  of  the  evacuation. 

Subroutine  INITIA 

INIT1A  is  called  by  subroutine  SITE.  The  purpose  of  this  subroutine  is  to  initialize  the  LLNL 
GRAFCORE  graphics  package  prior  to  the  plotting  of  the  road  network  and  associated  evacuation  in¬ 
formation.  INITIA  sets  up  the  reference  15-  by  15-mi  coordinate  system  used  for  the  plots. 

Subroutine  ISOPLT 

At  each  dispersion  update,  the  user  can  choose  to  plot  dose  isopleths  for  the  three  nuclides  1311, 133Xe, 
and  l37Cs.  The  contour  levels  are  supplied  by  Lawrence  Livermore  National  Laboratory  for  each  nuclide  at 
each  15-min  interval  and  are  stored  in  the  data  file  ACLEVS.DAT  on  unit  8.  These  isopleths  can  selec¬ 
tively  appear  on  the  screen  with  the  node  and  link  population  data,  which  makes  it  possible  to  determine 
the  parts  of  the  road  network  that  are  particularly  dangerous  as  far  as  exposure  is  concerned. 

Subroutine  LNKSET 

LNKSET  is  called  once  at  the  beginning  of  the  program.  It  determines  the  length  of  each  link,  the 
beginning  and  end  nodes,  the  grid  boxes  in  the  40  X  40  grid  through  which  each  link  passes,  and  the 
length  of  the  segment  of  the  link  within  each  grid  box.  The  length  and  location  of  the  segments  of  each 
link  are  required  for  the  dose  calculations. 

Subroutine  NODSET 

NODSET  is  called  once  at  the  beginning  of  the  program.  It  determines  the  grid  location  of  each  node. 

Subroutine  PLTEND 

PLTEND  is  called  to  terminate  each  plot,  both  road  network  plots  and  histogram  plots. 

Subroutine  PLTSET 

PLTSET  is  called  by  subroutine  ISOPLT.  At  each  dispersion  update  the  user  can  choose  to  superim¬ 
pose  current  dose  isopleths  on  the  plot  of  the  road  network.  However,  only  one  nuclide  can  be  displayed 
at  a  time.  PLTSET  asks  the  user  for  the  nuclide  of  interest,  loads  the  appropriate  dose  data  into  a  variable 
U,  dimensioned  40  X  40,  loads  the  appropriate  contour  levels  into  a  variable  ACLEVS,  dimensioned  4, 
and  returns  these  values  to  ISOPLT  to  be  plotted.  PLTSET  also  types  information  in  the  top  right  comer  of 
the  screen,  identifying  the  nuclide  of  interest  and  the  contour  levels  plotted. 


Subroutine  SITE 

SITE  is  called  at  the  start  of  the  calculation  and  at  each  dispersion  step  and  plots  site-specific  features 
as  a  prelude  to  the  plotting  of  isopleths  and  populations  designed  by  links  and  nodes.  These  features 
include  the  reactor  site  boundary,  the  reactor  building,  and  a  lake  within  the  site  boundary.  SITE  also 
plots  the  road  network  and  county  lines.  After  SITE  plots  the  links  and  nodes  of  the  road  network, 
subroutine  DISP  draws  over  the  links  and  nodes  to  indicate  population  density. 


Description  of  EVACD  Input  Files 

The  EVACD  program  uses  six  input  files  (Table  B-l)  as  described  in  detail  below. 


Table  B-l. 

EVACD  input  files. 

Unit 

File 

number 

name 

Contents 

8 

ACLEVS.DAT 

Contour  levels  to  be  plotted  each  time  step 

10 

EVACDB.INP 

Road  network  information 

11 

EVACDB.CAS 

Specific  evacuation  case  information 

15 

I131.LLNL 

40  X  40  doses  due  to  ,3,I  each  time  step 

16 

XE133.LLNL 

40  X  40  doses  due  to  l,3Xe  each  time  step 

17 

CS137.LLNL 

40  X  40  doses  due  to  u7Cs  each  time  step 

ACLEVS.DAT 

This  file  contains  contour  levels  to  be  plotted  for  each  nuclide.  There  are  18  dispersion  updates 
representing  15-min  intervals  and  four  contour  levels  for  each  nuclide  at  each  update.  On  each  line  in  the 
file,  the  first  four  numbers  are  contour  levels  for  ,31I  inhalation  dose,  the  next  four  are  contour  levels  for 
,33Xe  direct  dose,  and  the  last  four  are  countouT  levels  for  137Cs  inhalation  dose. 

EVACDB.INP 

This  file  contains  information  regarding  the  road  network,  plus  problem-specific  data.  The  first  line  is 
merely  a  title  describing  the  specific  case.  Line  2,  read  in  free  format,  is  as  follows: 

PPV,  VL,  T,  TPR1NT,  POPFRC,  NOPT1, 

where 


PPV 

VL 

T 

TPR1NT 

POPFRC 

NOPT1 


average  number  of  people  per  vehicle, 
vehicle  length  (mi), 
time  increment  (h), 

length  of  time  between  dispersion  updates  (h), 

fraction  of  population  which,  if  evacuated,  would  terminate  the  problem,  and 
an  option.  If  NOPT1  is  not  equal  to  1,  the  code  will  not  execute  subroutines  LNKSET  and 
NODSET  and  thus  will  not  set  up  the  geometry  necessary  to  calculate  population  doses. 
If  doses  are  desired,  NOPT1  should  be  set  to  1. 


Line  3,  also  read  in  free  format,  is 
NT1MES,  NNODES,  NLINKS,  NPOPN, 
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X 


where 


NTIMES 

NNODES 

SLINKS 

NPOPN 


maximum  number  of  time  intervals  that  the  problem  will  run, 
number  of  nodes  in  the  road  network, 
number  of  links  in  the  road  network,  and 

number  of  input  nodes,  i.e.,  nodes  at  which  people  can  enter  the  road  network. 


The  next  “NNODES"  lines,  read  in  format  15,  Al,  2F6.2,  are 


NODNM(N).  NODID(N),  XNODE(N),  YNODE(N),  N  =  l, NNODES, 

where 


NODNM 

NOD1D 

XNODE 

YNODE 


the  node  ID  number, 

the  node  type  ID,  with  "A"  denoting  an  input  node,  "Z"  denoting  an  output  node,  and 
'  B'  denoting  all  other  nodes, 
the  x-position  of  the  node,  and 
the  y-position  of  the  node. 


The  next  “NPOPN"  lines,  read  in  free  format,  are 


ID,  POPINUD ),  CAPNODOD),  PIN(ID), 


where 


ID 

POPIN 

CAPNOD 

PIN 


the  ID  number  of  the  input  node  and  is  equivalent  to  NODNM  above, 
the  initial  number  of  people  at  the  input  node, 

the  capacity  of  the  input  node,  i.e.,  the  number  of  vehicles/h  that  can  enter  the  input 
node,  and 

preference  factor  (fraction  from  0-1)  for  the  input  node. 


The  last  " NLINKS "  lines,  read  in  free  format,  are 


LID(L),  LLD(L),  LNUMUL),  FFS(L),  CAP(L),  AL(L),  ARID,  PFIN(L),  GS1N(L).  LNIN(L),  LNOUT(L), 
QINIT(L),  VEHINUL),  L=  1,  NLINKS. 

where 


LID 

LED 

LNUML 

FFS 

CAP 

AL 

AR 

PFIN 


CSIN 

LNIN 

LNOUT 

QUNIT 

VEHINI 


the  link  ID  number, 

the  length  of  the  link  (mi), 

the  number  of  traffic  lanes  for  the  link, 

the  free  flow  speed  of  the  link  (mi/h), 

the  capacity  of  the  link  (vehicles/h);  this  is  the  maximum  number  of  vehicles  that  have  a 

reasonable  expectation  of  passing  over  a  given  section  of  roadway  in  1  h, 

the  left  turn  adjustment  factor  for  the  link  (currently  not  in  use), 

the  right  turn  adjustment  factor  for  the  link  (currently  not  in  use), 

the  preference  factor  (fraction  from  0  -  1)  for  the  link  into  its  output  node,  i.e.,  the 

fraction  of  the  time  that  traffic  from  the  link  can  pass  through  the  output  node,  relative 

to  other  links  with  the  same  output  node, 

the  green  split,  or  fraction  of  the  time  that  traffic  signals  on  the  link  are  green;  if  a 

negative  number  is  used  the  program  will  use  the  green-split  variable, 

the  ID  number  of  the  input  node  for  the  link, 

the  ID  number  of  the  output  node  for  the  link, 

the  number  of  vehicles  initially  in  a  queue  at  the  end  of  the  link,  and 

the  initial  number  of  vehicles  on  the  link. 


A 


EVACDB.CAS 

This  file  contains  information  regarding  specific  parameters  for  the  evacuation  case  being  run.  The 
first  line  is  merely  a  title  describing  the  case.  Line  2,  read  in  free  format,  is  as  follows: 

TNOTIC,  TDELAY,  TUPD, 

where 

TNOTIC  =  the  notification  time,  or  the  time  that  elapses  following  the  release  until  people  are 
notified  to  evacuate  (h), 

TDELAY  —  the  time  that  elapses  following  notification  before  people  begin  to  evacuate  (h).  and 
TUPD  =  the  time  between  dispersion  updates  (h). 

Line  3,  read  in  free  format,  is 

XZERO,  YZERO,  NR1NGS,  NZONES,  . 

where 


XZERO  =  the  x  position  of  the  release  point, 

YZERO  =  the  y-position  of  the  release  point, 

NR1NGS  =  the  number  of  "rings"  around  the  release  point,  established  for  the  purpose  of  watch¬ 
ing  the  progress  of  the  evacuation;  in  the  sample  input,  ten  1-mi  rings  were  used,  with 
an  eleventh  ring  from  10  mi  out  to  15  mi  for  the  evacuated  population,  and 
NZONES  =  the  number  of  evacuation  zones  under  consideration. 

Line  4,  read  in  free  format,  is 

R1NG(N),  N  =  1,  NR1NGS. 

where 

RJNG(N)  —  the  outer  radius  of  the  nth  ring  (mi). 

Following  line  4,  the  file  contains  "NZONES"  pairs  of  lines,  as  follows: 

ZONE(f),  NSEC(I) 

MSECS(J.l),  /  =  1, NSEC(J), 

where 


ZONE(f)  =  the  outer  radius  of  the  ;th  evacuation  zone  (mi), 

NSEC(J)  —  the  number  of  sectors  included  in  the  ith  evacuation  zone,  and 
MSECS(J.l)  —  the  ID  numbers  of  the  sectors  included  in  the  ;th  evacuation  zone. 

Usually  the  area  around  the  release  point  will  be  divided  into  16  equal  sectors  numbered  1  through  16.  It  is 
possible  to  set  up  an  evacuation  area  in  a  shape  other  than  a  circle  by  having  an  outer  evacuation  zone 
consist  of  different  sectors  from  an  inner  zone.  For  example,  a  key  shaped  area  can  be  established  by 
having  an  inner  zone  out  to  5  mi  consist  of  all  16  sectors  and  have  an  outer  zone  out  to  10  mi  consist  of 
about  3  sectors.  The  sample  problem  uses  only  one  zone. 

The  last  line  of  the  file,  read  in  free  format,  is 

NGR1DX,  NGR1DY,  XCORN,  YCORN,  DELX,  DELY, 

45 


where 


NGR1DX 

NGRIDY 

XCORN 

YCORN 

DELX 

DELY 


the  number  of  grid  boxes  in  the  x-direction  for  which  dose  information  is  given, 
the  number  of  grid  boxes  in  the  y-direction  for  which  dose  information  is  given, 
the  x-position  of  the  lower  left  comer  of  the  grid  (mi  from  the  center), 
the  y-position  of  the  lower  left  comer  of  grid  (mi  from  the  center), 
the  length  of  a  grid  box  in  the  x-direction  (mi),  and 
the  length  of  a  grid  box  in  the  y-direction  (mi). 


I131.LLNL,  XE133.LLNL,  CS137.LLNL 

These  three  files  contain  the  dose  information  for  the  nuclides  u,l,  133Xe,  and  137Cs,  respectively.  They 
are  set  up  on  the  basis  of  a  41  X  41  grid,  with  41  X  41  values  at  each  of  18  times.  The  EVACD  program 
uses  a  40  X  40  grid,  so  it  ignores  data  corresponding  to  the  41st  row  and  41st  column.  The  order  of  the 
data  is  as  follows: 


7=1, 

/=1  +  41' 

J=2, 

J=1  +  41 

’ 

7=41, 

/=  1  +  41. 

7=1, 

/=  1  -f  41' 

7=2, 

f=l  -f  41 

► 

7=41, 

f=l  +  41. 

7=1, 

/=1  +  41' 

7=2, 

/=*1  +  41 

. 

7=41, 

/=1  +  41,. 

Time  since  release  =  0.25  h  , 


Time  since  release  =  0.50  h  , 


Time  since  release  =  4.50  h  , 


where  /  refers  to  the  y-coordinate  of  a  data  point  and  I  refers  to  the  x-coordinate  of  a  data  point.  Data  are 
read  in  12E10.2  format. 
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Appendix  C: 

Listing  of  EVACD  Main  Routine, 
Subroutines,  and  Common  Files 


program  evacdb 
real-time  evacuation  model 

FEBRUARY  1901 


INCLUDE  'DSKD: PARAM.COM' 

INCLUDE  ’ DSKD : NODE . COM ‘ 

INCLUDE  ’DSKD: LI NK.COM' 

INCLUDE  -DSKD: DOSE. COM' 

INCLUDE  'DSKD -.ROAD. COM’ 

Include  " dskd : t ex t . com  * 

DIMENSION  T I TLE  <  20 ) . ET I TLE (20) 

DIMENSION  NODNM ( N 1 > .POPIN(NI ) .CAPNOD(Nl ) . 

*  P I N »  N l ) . NLFNlNl > ,NLTO<Nl  ■ ,VEH1N(N1 ) ,VEH(Nl ) . 

*  VEHND(N1  l.MNl  »,I’F(NI  )  .NODOUTiN!  )  ,EXlT(Nl  ) 

DIMENSION  LID(N2> .LLDIN2) ,FFS(N2) . DJAM ( N2 ) . 

*  PF1N(N2) .GS1N(N2> ,Q1N1T(N2> . 

*  VEHIN1  ( N2  )  ,  VEH I CS (  N2  >  .QUE0UT(N2)  .QUELFN(N2> . 

*  Ml NP  C  N2 ) .DENSIN2) ,SPEED(N2> .FLOV(N2> .  V EHTN ( N2 > , 

*  VLEFT(N2> -EXVEH(N2> .CS'N2) . CAPLNK  (  N2 ) .CAPiN2> . 

*  AL<  N2 ) .  AR  <  N2 ) ,APPCAP(N2) . V  I t N2 ) . P( N2 ) . PT( N2 > . 

*  VO<  N2  » . MOUT ( N2 ) . W <  N2 > . P  M  N2 ) 

DIMENSION  LINKIMN1  . N3 ) . LI NKOT( N 1  .N3) 

RE4L  LLD . LENCTH . NODOUT . M I NP , MOUT 
DATA  A.Z/'A  ' . 'Z 
KVARL* l2*N2+a*Nl*2 

open  (  uni  1=8,  slat  us  *  *  o  I  d  ’  .  f 1 1 e = ‘ dsk : ac 1 fvj . da t ’  ) 

OPEN (UNIT* I®. status* ’old’  .file* ’ ds>. : EVACDB . I NP -  > 

OPEN (UN  IT* 11  .status* 'old ' .flic* ' dak : EVACDB . CAS '  ) 

OPEN  ( UN  IT*12. status*  ' uev, '  ,  f  1 1  e *  '  dsk  :  EV  VCDB  .OUT '  > 

OPEN (UN IT* 13. status* 'old ' .file* '  tsk s 1131 . LLL ' ) 

OPEN ( UN  IT* 1 6. s ta tus* 'old ’ , f i I e= ' dsk : xc ! 33 . LLL ' > 

OPEN (VNIT=17. status* 'old ’ . f i le= ' dsk : cs 1 37 . LLL ' > 
c  TYPE  3000 

3000  FORMAT (  I H I  ,  *  BE;C IN  EVACD  EVACUATION  PROCRAM ’  . 

*  ✓✓  '  APPROPRIATE  FILES  ARE  BEING  READ  IN') 

READ( 10. 1  1 06  IT I TLE 

1106  FORMAT (20A4) 

REAP (10.*) PPV . VL .T . TPR I  NT . POPFRC . NOPT 1 
REAP(  lO.*)NTIMES,NNODES.  NLI NKS  . NPOI’N 

READdO.  1  105  )  ( NODNM ( N  )  . NODI  D(  N >  . XNODEi N)  .  YNODE<  N )  , N*  I  ,NNODES> 
c  TYPE  1300. TITLE 

VR1TE( 12.1 600 1TITLE 

1600  FORMAT  < 20A4 ) 

150©  FORMAT! IHI , 10X.2QA4) 

VRITE< 12.1601 ) PPV , VL ■ T . NTl MES . TPR I  NT . NNODES .NLI NKS . POPFRC . NOPT  1 

1601  FORMAT' 3E1 0.2. 1 6. El  0.2. 2 1  6. El  0.2 , lb) 

c  TYPE  1501 .PPV. VL.T, NT  I MES, TPR I  NT. NNODES. NLI NKS. POPFRC. NOPT 1 

1501  FORMAT (//'2X,  *  INPUT  PARAMETERS'  . 


* 

z'/'SX ,  ’  PEOPLE  PER  VEHICLE 

=  * .3X.F10. 1  . 

* 

✓3X. 'VEHICLE  LENCTH' Ml) 

=  '  ,3X. IPE10.2. 

* 

✓5X .’TIME  I NCREMENT ( HR  * 

=’ , 3X . OPF I0.3. 

* 

/5X. ’NUMBER  OF  TIME  INCREMENTS 

= ’ ,7X, 16. 

* 

/5X. 'REPORTING  TIME(HR) 

= ’ .3X.F10.3, 

* 

/5X. ’NUMBER  OF  NODES 

=' .7N.16. 

* 

/5X, ’NUMBER  OF  LINKS 

*  ’  . 7X . I 6 . 

* 

✓5X . ’CUTOFF  FRACTION 

*’  .3X.F10.4 

* 

/'5X,  ’NOPTI 

*  •  ,7X, 16) 

tt03  FORMAT* 15, Al .2E6.2) 

DO  1  N=l .NPOPN 

READ  <  1 0  .  * )  I D .  POP  1 14  C  ID)  .CAPNOD(ID)  .PINUD) 
C  TOTAL  UP  ALL  POPULATION  ENTER  I NC  'A'  NODES 
TTPOPN*TTPOPN»POPI N  < ID) 


Figure  C-l.  EVACD  main  routine. 
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I  CONTINUE 

READ! |0.*HL1D(L) . LLD(L> . LNVML(L) .  FFS(L> .CAP(L) . AL( L ) , AR ( L ) , 
*  PFIN<L>  .CSIN(L)  . LN 1 N i L  >  .LiOUTlL)  .OIN1TIL) . 


*  VEH1NI  <L)  .  L  -  1  . NLINKS) 

C  READ  IN  EVACUATION-SPEC  IF  1C  DATA 

READ! I  I . I  106  IETITLE 
READ( 1 1 ,*)TNOTIC .TDELAY ,TUPD 
READ(1  I  , * )XZ£RO ,  YZERO .  NR  1  NCS  ,  NZONES 
READ( 1  I .*)<R1NC(N> ,N=I .NH1NCS) 
e  TYPE  3001 

3001  FORMAT < 1R1 , '  SPECIFIC  EVACUATION  CASE  IS:*) 

TYPE  2999 ,  ET1TLF. 

2999  FORMAT! /1X.40A4) 

c  TYPE  3002 

3002  F0RMATi//5X. 'SPECIFIC  PARAMETERS  FOLLOW:’) 

c  TYPE  3003 .TNOTIC .TDELAY .TUPD 

3003  FORMAT*/ I  OX. 'NOTIFICATION  TIME (HR)  =’.F6.2. 

*  /IOX.  •PREPARATION  TIMFOIK  >  =  ’.F6.2, 

*  /I OX. ’DISPERSION  UPDAT£(K!l>  =  ’.F6.2) 

c  TYPE  3004 . NZONES 

3004  FORMAT! /\ OX. ’NUMBER  OF  EVACUATION  ZONES  *’.12) 

DO  200  J*l . NZONES 

READ! 11 . * )ZONE! J ) <  NSEC  1 J ) 

NS=NSEC(J) 

READ! 1 1 . * ) ! MSECS ( J . I) . I  -  I .NS) 
c  TYPE  3005 . J , ZONE! J  > . ! MSECS  <J. I). 1  =  1. NS) 

3005  FORMAT! //I OX. ’ZONE  NUMBFR’.12.’  IS  OUT  TO ’ ,F5 . 1 . ’  MILES’. 

*  / I 0X, ’ INCLUDING  SECTORS: 

*  X10X.I613) 

200  CONTINUE 

READ Ill.*) NCR  I PX . NCR  I DY . XCORN . YCORN . DELX . DELY 
IF1N0PTI  .NE.  1 ) GOTO  05 

C  CALL  NODSET  TO  SET  UP  NODE  LOCATION i RELATI VE ) 

CALL  NODSET 

C  !F  ERROR  DETECTED  COTO  END 

I F !  NFFiHOR  .  F.0  .  I  )CO TO  171 

C  CALL  LNKSr'.T  TO  SET  UP  LINE  LENGTH  AND  LOCATION  DATA 
CALL  LNFSET 
03  CONTINUE 

C  TOTAL  VP  ALL  POPULATION  TO  BE  EVACUATED 
DO  201  Nl I . N NODES 
I F ! NODI  DIN)  .NE.  ’A  ')COTO201 
I T I NODVAC ! N )  NE.  1  .AND.  NOPTl  .  EQ .  I ) COTO  201 
TPOPN =TPOPN +POP ININ) 

20 1  CONTINUE 

I F< TPOPN  . NE.  O.ICOTO  172 
TYPE  2000 

2000  FORMAT! ’  PROGRAM  TERMINATED  FOR  ZERO  POPULATION  EVACUATION’) 
COTO  171 
172  CONTINUE 

C  CALL  CAPACY 

C  SET  UP  NODE  AND  LINK  INFO 
DO  120  L«l. KLINES 
C  CALCULATE  JAM  DENSITY 

DJAMIL) *4 .*CAP(L)/FFS!L> 

C  TOTAL  UP  ALL  POPULATION  ON  LINKS 
TPOPL=TPOPL+VEH 1 N I <L)*PrV 
NOUT^LNINIL) 

NIN=LNOUTIL) 

NLFN ININ) - NLFN I N I N ) + 1 
NLFO ( NOUT ) * NLFO I NOUT  ) ♦ 1 
LI =NLFN ININ) 

L2= NLFO (NOUT) 

LINKIN'NIN.Ll )-LlD(L) 

LINK0T(N0UT.L2)=LID<L) 

120  CONTINUE 

1100  FORMAT (2X, I5.4X.A1 ) 

1101  FORMAT! 5X , ’ IN  =  ’.513) 

1102  FORMAT I 3X . ’OUT  «  *,515) 
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c  CALL  LNKPLT 

C  TOTAL  UP  POPULATION  INITIALLY  AT  *A*  NODES  AND  ON  ALL  LINKS 
TOTPOP»TPOPL+TPOPN 

c  TYPE  I6O4.T0TP0P.TP0PN.TO1TL 

1604  FORMAT (//'  TOTAL  POPULATION  FOR  EVACUATION  **.F10.0, 

*  ✓  *  TOTAL  POPULATION  AT  NODES  =*.F10.0. 

*  /’  TOTAL  POPULATION  ON  LINKS  =’,Fie.O) 

C  CALCULATE  INITIAL  NUMBER  OF  VEHICLES  VAITINC  TO  ENTER  INPUT  NODE 
DO  5  N=l , NNODES 
VEH 1 N ( N ) =  POP  I N ( N ) /PPV 
C  VR I TE (12.1 103 ) N .VEH 1 N<  N) .POPIN(N) 

5  CONTINUE 

1103  FORMAT ( IX. 'FOR  NODE’. 15, •  INITIAL  VEHICLES  =  *.F6.0. 

*  •  INITIAL  POPULATION  =  '.F6.0I 

C  TOTAL  TIME  DELAY  BEFORE  EVACUATION  START  IS  TNOTIC+TDELAY 
TOTDE=TNOTIC+TDELAY 
c  TYPE  3006 

3006  FORMAT ( ’  BECIN  RELEASE  EVENT  AT  TIME  =  O’) 

THDATA=0. 

T I  MOO 

C  START  TIME  INCREMENT  LOOP 
ndone -  0 
TNECT-TUPD 
DO  10  HT*1  .  NTIMES 
ndoser=0 
NTM=  NT 

I F (TMPATA  .LT.  -.001  .OR.  TMDATA  .CT.  .OODCOTO  71 
nJos«r: 1 
NTM1N-NT-I 
TYPE  3007. TIME 

3007  FORMAT <1 HO. 5X. ’TIME  SINCE  RELEASE! HR)  -’.F6.2) 

ISITRD=NT-1 

75  CALL  SITE! ISITRD. TIME) 

I F ( NT  .EQ.  1 )COTO  72 
CALL  DISP(NTM.VL) 

CALL  I SOPLT(TNEXT) 

CALL  l’t.TEND 

73  TYPE  9000 

9000  FORMAT! ✓ *  DO  YOU  VANT  TO  SEE  ISOPLETHS  FOR  ANOTHER  NUCLIDE* 

I  .*  AT  THIS  TIME  STEP?  V '  ENTEH  FOR  YES  OR  *N"  FOR  NO.*) 

ACCEPT  9001. IANS 

9001  FORMAT! Al) 

Irln  *  a  1 r*upca«e ( Udesor ( 1 ans ) ,%dc scr ( 1  ins ) ) 

IF( 1ANS.EQ. *Y* )  GOTO  75 
1F< IANS.EQ. *N* )  GOTO  74 
CALL  ERROR 
GOTO  73 

74  CALL  HISTO(RINPOP. time) 

•next1 tneit+tupd 

72  CONTINUE 

CALL  PLTEND 

I F  <  NDONF. .  F.Q .  1  )  GO  TO  11 
CALI.  DATA 
TMDATA = -TUPD 
71  CONTINUE 

TMD AT  A  r  TMDAT  A+T 

C  CHECK  TO  SEE  IF  EVACUATION  HAS  STARTED 
T I  MF.»  FLOAT  ! NT ) »T 
C  TYPE  * ,TI ME 

IFITIME  .CE.  TOTDE)NTGO=NTCOeI 
TOTEX-O . 

C  START  LINE  LOOP 

DO  20  L= I .NLINKS 
IFl NT  .NE.  I )COTO  21 
C  SET  UP  INITIAL  LINK  OIT  QUE 
VEH I CS ( L )  =  VEH I N I  (L) 

QUEOUT <L)*QINIT!L) 

COTO  19 

C  CALCULATE  LENCTH  OF  OUT  QUE 
21  CONTINUE 

VEHICS(L)«VLEFT<L)+MINP<L) 
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19  QUELEN ' L  )  *QULOUT ( L) *VL/LNUML( L ) 

c  TYPE  1 203 , NT . L . VEH 1  CS ( L  > 

LENCTH  =  LLD  ( L )  -Ql'ELEN  ( L ) 

I F ( LENGTH  .CT.  0)COTO  22 
LENGTH -0 . 

DENS(L>=0. 

COTO  20 

22  CONTINUE 

C  CALCULATE  DENSITY! PER  LANE) 

DENS  (  L  )  =  VEH  I CS  ( L )  /LENCTH/LNVML  ( L ) 

23  CONTINUE 

C  CALCULATE  FLOW  SPEED 

STEED'  L)  =  FFS(L)*'<  1  . -DENS ( L>/DJAM ( L>  ) 

C  CALCULATE  LINK  FLOW 

FLOV(L'=DENS(L>*SPEED<L)*LNUML(L> 

C  CALCULATE  VEHICLES  REACHING  END  OF  LINK  OR  OUT  QUE 
VEHTN < L ) =FLOW ( L  >  *T 

I F ( VEHTN ( L)  .CE.  VEH I CS (L) ) VEHTN < L > = VEH I CS (L) 

C  VEHICLES  REMAINING  ON  LINK 

VLEFT <  L ) =  VEH I CS  <  L ) -VEHTN ( L ) 

XL(L. 1 )  =  VEH I CS  < L) 

XL! L .2  >  =QUEOUT ( L ) 

XL! L ,  3  >  =  LENCTH 
XL(L,4>-DENS<L) 

XL(L.5'=SPEEU(L» 

XL!  L .  O  >  =  FLOW  (  L  1 
XL ( L . 7  »  =  VEHTN ( L » 

XL ( L . 8 1 =  VLEFT  !  L  > 

GUF.OUT '  L )  sQUFOUT!  L)  ♦ VEHTN  ( L ) 

QUELEN  •  L  >  *QUEOUT<  L)  *VL-'LNUNL(  L) 

LENCTH  =  LLD ( L ) -01  ELEN ( L ) 

XU  L .  9  >  iQUEOUT (  L) 

XL ( L .  1 0  >  =  LENGTH 

IF! LENCTH  .LE.  0>LENCTH=0. 

DENS ! L ) = VLEFT ! L i /LENCTH^  LNUML ! L > 

EX VEH ( L  >  =  LENGTH* ( DJ AH ! L » -DENS ! L  > >*LNUML(  l. ) 

I F !  F.XVEH  ( L )  .LE.  0  .  )EXVEH(L)  =0  . 

C  WHITE! 12. I  1 04 ) NT . L . VEH I CS ! L ) .QUELEN(L) .DENSIL) ,SPEED!L) .FLOWIL) . 

C  *  VEHTN(L) ,VLEFT(L' ,£XVEH(L> 

1I©4  FORMAT' IX. 215. 8F7.0) 

20  CONTINUE 

END  OF  LINK  LOOP 
zzin=0.  I pop  entering 

zzout=P.  -pop  leaving 

zzleft=0.  Ipop  on  links 

zzque-O.  ! queue  on  links 

do  8070  l=l.nlinks 
zzqiie  -  zzque  +  xl  (  1  .  9  ) 

8070  zz lef t  =  zz left *x 1(1.8) 

do  8071  n=l.nnodes 
zz I n  =  zz I  n+m  (  n  .2  ) 

6071  zzon t * zzou t +xn ( n , 3 ) 

zz  t  o  I  =  zz 1 e  f  t  -*-zzque  +  zz  1  n  +  zzou  t 
call  dscalc(nta.t.tupd. ndoser ,ppv) 

C  START  NODE  LOOP 

DO  30  N  =  I  , NNODES 
C. CHECK  FOR  IN  OR  OUT  NODE 
N0DT=0 

IF! NODI D( N )  .NE.  ’A  ’ >GOTO  31 
IF(NODVAC(N)  .NE.  1  .AND.  NOPTl  .EQ.  I ) COTO  31 
NODT  = I 

C  SET  UP  INPUT  INTO  NODE 

I F ( NTGO  .EQ.  1 >VEH(N)*VEHIN(N> 

C  MAXIMUM  INPUT  IN  VEHICLES  PER  TIME  INTERVAL 
YEHND ( N ) »  CAPNOD  <  N ) *T 

I F ( VEH ( N )  .LE.  VEHND(N) )VEHND(N)*VEH(N) 

C  TYPE  1 200, NT. N. CAPNOD! N) .VEH(N) .VEHND! N) 

1200  FORMAT! IX. 215. 3F7.0) 

31  CONTINUE 

XN(N, 1 )*VEH(N> 

C  TYPE  1201 .NT.N.NODID(N) 

1201  FORMAT! IX. 213, 5X.A1 ) 

I F ! NODI D(N )  .EQ.  ’Z  * )N0DT*2 
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C  SET  UP  NUMBER  OF  INPUT  AND  OUTPUT  LINKS  FOR  NODE 
NL1 ®NLFN(N) 

NLO=NLFO(N) 

V(N)*0. 

C  LOOP  OVER  INPUT  LINKS 

IF ( NLI  .EQ.  0 )GOTO  46 
DO  32  L: I .NLI 

C  LINK  IS  INPUT  LINK  L  TO  NODE  N 
LINKELINK1N(N.L) 

VW  ( LI  Ni: )  *QUEOUT(  LI  NK ) 

V ( N  >  =  W ( N  > +VW (LINK) /LNUNL (LINK) 

32  CONTINUE 

C  ADD  POTENTIAL  INPUT  VEHICLES  FROM  INPUT  NODEITVPE  A) 
I F ( NODT  .EQ.  1 )V(N)=VMN)+VEHND<N> 

IF(MN)  .LE.  0.  )V(N)  =  I  . 

C  CALCULATE  GREEN  SPLIT 
DO  33  L= I .NLI 
LIN»LINKIN(N.L» 

I F  ( GS I  N  ( L I  NK  )  .C-E.  0.)COTO34 

cs«lini:)=W(link)/lnunl'Link)/v(n> 

GOTO  35 

34  CS<  LI NK ) =CSIN(LINK ) 

35  CONTINUE 

C  CALCULATE  ENTRY  INPUT  GREEN  SPLIT 

CAPLNTILINF.)  =CAF<LI  NK  >*AL<LI  NK  )*AR«  LINK) 

C  USE  FIRST  CONSTRAINT 

APPCAP  -LINK)  =CAFI.NK  (LINK)  *GS  (LINK) 

VI (LINK ) =APPCAP( LI NK )*T 

IF(V1  (LINK)  .CE.  VV(LINK))VMLINK)*VV(LINK) 


33 

CONTINUE 

I F ( NODT  .EQ. 
GOTO  47 

1 )ENTGS=VEBND(N)/V(N) 

46 

CONTINUE 

I F < NODT  .EQ. 

1 )ENTCS= 1 . 

47 

CONTINUE 

C  CALCULATE  OUTBOUND  LINK  PREFERENCES 
PF  f  N> -0. 

IF(NODT  .EQ.  2 )GOTO  41 
DO  36  L= 1 .NLO 
LINI=LINKOT(N.L) 

PF( N ) =PF(N ) +PFI NIL INK )*SPEED( LINK) 

36  CONTINUE 

1 F ( PK ( N )  -LF..  0 .  )PF ( N )  =  1  . 

DO  37  L= 1 . NLO 
LIN»LINKOT(N.L) 

P«LINK)  =  PFIN(LINK)*SPEED(LINK)/,Pr(N) 

37  CONTINUE 

DO  38  L 1 = I , NLO 
LI NVI *LINKOT(N ,LI ) 

PT(LINK1 ) =0 . 

I F ( NLI  .EQ.  OIGOTO  48 
DO  39  L2= I .NLI 
LINK2=LINKIN(N.L2> 

P I ( L I NK2 ) 1 1 . 

PTILINKl )»PT<LINKl )*VI <LINK2)*P(LINK1 J*PI<LINK2> 

3*  CONTINUE 

IF (NODT  .EQ.  l)PT(LINKl)«PT(LINKt)^VEBND<N)*P(LINKt)*PIN(N> 
GOTO  49 

48  CONTINUE 

I F I NODT  .EQ.  I )PT(LINK1)«VEHND(N) 

49  CONTINUE 

C  CALCULATE  OUTBOUND  LINK1  VOLUME  RECEIVED (SECOND  CONSTRAINT) 
VOCLINKl >«CAP(LINK1 )*T 

1FIPT ( LI NK1 )  .LE.  V0(LINK1 ) )VO(LINKI )*PT(LINK! ) 
1F<EXVEH(LINK1 )  .LE.  VO(LINKI > )V0(LINKI >*EXVEH(LlltKI ) 
MlNPILINXl >«0. 

1F(PT(LINKI )  .LE.  • . )PT(LINKI )■ 1 . 

IF( NLI  .EQ.  9>COTO  42 
DO  40  L2«l .NLI 
LIN>:2EL!NKIN(N.L2) 

PI ( LINK2 ) *  I . 

H0UT(LINK2)«V1 (LINK2)«P(LINK1 )*PI (LINK2)*V0(LINKI )/PT<LINKl ) 
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MI  NP ( LI NK 1 >«MINP(LINKl )+N0UTtLINK2) 

43  CONTI NUE 

QUEOUT  v  L I NK2  )  *QUEOUT  ( L I NK2 )  -MOUT  ( L 1 NX2 ) 

I F (QUEOUT ( LI NK2  *  .LE.  O . >QUEOUT ( LI NK2 > *0 . 

XL<  LI  NK2 . 1  1  >*N0VT(LINK2  ' 

XL<LINK2. 12)=QUE0UT(LINK2) 

40  CONTINUE 

42  CONTINUE 

in  HOOT  .NE.  OC0T0  44 

NODOUT  N)=YEHND<N)*P(LINY.l  )*? IN(N)*VO(LINKl  >/PT(LINKt  > 
MINPd.INKI  >*NINP(LINKl  >-NODOUT(N) 

VEK <  N ) = VEH l N ) -NODOUT  ( N ) 

IF(VEH'N)  .LE.  0. >VEH(N)*0. 

44  CONTINUE 

C  TYPE  1202. NT. N. LINK! ,MIKP(LINK1 ) .NOOOUTIN) .VEH'N) 

1202  FORMAT- IOX. 315, 3F7.0) 

38  CONTINUE 

COTO  52 

41  CONTINUE 

DO  43  L* 1 .NLI 
L1NE*LINK1N(N.L> 

EXIT<N  *  =  EXIT(N)*V! (LINK) 

C  TYPE  1 203 .NT , N . EXIT*  N ) 

1203  FORMAT! I0X.215.F7.1 ) 

MOUT ( LI NK  >  *VI (LINK) 

QUEOUT  (LINK)  *OUl.OUT  (-INK)  -MOUT  (LINK) 

XL < L 1  NK . 1 1 >*HOin  (LINK) 

XL(L1T;K.  12>*QUE0UT(LJNK> 

45  CONTINUE 

52  CONTINUE 

C  I E ( NODT  .EQ.  DTYPE  1203  .NT . N .VEH (N) 

I F < NOW  .EQ.  1  )XN(N.2)«VEH<N) 

I F ( MOOT  .EQ.  2 )XN (N.3)=EX1T(N) 

TOTEX  * TOTEX +EX I T ( N  > 

30  CONTI  NUF. 

end  OF  NODE  LOOP 
FRC*TOTEX^TOrrOP*PPV 
c  TYPE  1  605 .  NT .  FRC 

1603  FORMAT  ’  FOR  INTERVAL’ .  14.  ’  FRACTION  OF  EVACUATED  POP  «’,F6.4 
C  VR1TE(2C'NT)XL.XN. TOTEX. FRC 

I F ( NTM  .EQ.  I )CALL  D I SP  ( NTM . VL ) 

I F <  FRC  .LT.  I’OPFRC  )GOTO  10 

LTINE=NT 

TYPE  1405.LTIME 

1405  FORMAT  <  5X , ' LAST  TIME  INTERVAL  IS  NUMBER  ’,15) 

ndone  * 1 
COTO  73 

10  CONTINUE 

END  OF  TIME  LOOP 
LTI ME=NTIHES 

11  CONTINUE 
CALL  JECCXX 
WRITE) 12,1 602 )LT! ME 

1602  FORMAT' 16) 

WRITE) 12.1612) (NODID(N ) ,N* 1 .NNODES' 

1612  FORMAT -80 A I ) 

TYPE  1603, LTIME, FRC 

1603  FORMAT)  ✓v’  LAST  TIME  INTERVAL  *’.16. 

*  /'  FRACTION  OF  POPOUT  *’.F6.4) 

171  CONTINUE 

END 
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SUBROUTINE  DATA 

C  SUBROUTINE  TO  READ  APPROPRIATE  DOSE  DATA 
C  PRESENTLY  THREE  FILES  EXIT 
C 

c 

INCLUDE  ’DSHD s DOSE. COM' 

DIMENSION  A(41 ) 

DO  20  J=l .4© 

REAIM  1 5 ,9001 )  ( DOSEI I . J . 1 > . 1 = 1 ,40 ) 
READ! 16,9001 )  ( DOSE ( I . J .2) . I = 1 .40 > 

20  READ( 17.9001 )  (DOSE< I . J .3) . I = 1 ,40 ) 
READ! 15.9001 )  A 
READ! 16.9001 )  A 
READ* 17.9001 >  A 
9001  FORMAT (1 2E10 . 2  > 

RETURN 

END 


Figure  C-2.  Subroutine  DATA. 


SUBROUTINE  DISP(NTM.VL) 

C  SUBROUTINE  TO  DISPLAY  LINK  AND  NODE  CURRENT  NUMERICAL  DATA 

C 

C 

INCLUDE  'SYSOLIBRARYsCOLORDEF' 

INCLUDE  'DSKDsPARAM.COM' 

INCLUDE  'DSKDsNODE.COM' 

INCLUDE  'DSKDsLlNK.COM' 

INCLUDE  'DSKDsROAD.COM' 

DIMENSION  XLMONE ( N2 . 1 2 ) .XNMONEINl ,3) ,X(2> ,Y(2) 

DATA  DENMAX . W I DMAX/500 ..1.2/ 

I ROAD* YELLOW 
IQUE’ORANCE 
N COUNT* 0 

JFINTM  .NE.  1 JCOTO  30 
I  CONTINUE 

DO  10  L= 1 , NLINKS 
DO  10  J*1  .  12 

10  XLMONF.<L.J)*XL<L.J) 

DO  20  N* 1 . NNODES 
DO  20  J* I ,3 

20  XNMONE ' N , J ) *XN ( N . J ) 

IFINCOUNT  .EQ.  DCOTO  100 
30  CONTINUE 

DENCOP = V I DMAX/DENMAX 

3000  FORMAT < 1H 1 , •  LINK  DENSITIES  ARE  AS  FOLLOWS:’. 

*  //12X, ’DF.NSITY(VEH/MI ) '  . 

*  /2X. ’LINK’ ,4X. ’FRESENT* .6X. 'PAST' ) 

DO  40  L=1 .NLINKS 

DENS I T* XL < L . 1 2 ) »LNUML < L 1 
e  !F(L  . EO .  50  .OP.  L  .EQ.  I001PAUSE  21 

V I DTH  =  DENS I T*DENC0R*2 . 0 
IF1WIDTII  .CT.  WJDMAX)WlDTH=WTDrtAX 
CALL  JSLNSZI WIDTH) 

F1N=LN1N<L) 

NOUT*  LHOUT  ( L ) 

Xi 1 »*XNODE(NIN) 

X12i*XN0DE(N0UT> 

Yi 1 )*YNODE<NIN> 

YI2)*YN0DE(N0UT) 
c  denelty  color** 

CALL  JSCRXXI I ROAD) 

IRTN* JPPL2AIX.Y.2) 

XT*X( I )-X(2) 

YT*Y( 1 >-Y<2) 

RLEN'SORT ( XT*XT*YT*YT ) 

OLEN*XL(L. I2)*VL/LNUML(L> 

RATIO*OLEN/RLEN 


Figure  C-3.  Subroutine  DISP. 
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XQ=X<2 )+XT*RATIO 
YQ=Y<2>*YT*RATIO 
CALL  JSCRXX(IQUE) 

X< 1 •  *XQ 
Y<  I  >  *  YO 

OW|  R* FLOAT? LNUML(L)  )*V1DMAX 
CALL  JFLNSZlGVIIU 
IRTN  =  .;|-PL2A(X.Y.2) 

40  CONTINUE 

3001  FORMAT ( I  5 ,4X , F7 .2 .4X ,F7 .2 ) 

c  PAUSE  22 

3002  FORMAT? I  HI . *  NODE  INFORMATION  FOLLOWS:', 

*  /V13X .  'Ql'E<  VEHICLES)  •  ,4X.  'EXITINCfVEHICLES)  *  . 

*  ^2X . 1  NODE • . 4X . ’ PRESENT • . 6X . ’ PAST  * . 

*  4X , *  PRESENT  ’ , 6X , ' PAST ' ) 

CALL  JSLNSZ'O) 

DO  50  N  = 1 .NNODES 

1F<N0DID(N>  .EQ.  'B  * )COTO  50 

IF(N0DID<5)  .EO.  ’A  ’ )COTO  51 

IF<  NODID(N)  .EQ.  ’Z  ’» GOTO  52 

51  IF<XN<N.2>  .LE.O.O)  GOTO  50 
XLOC*ALOC 10(XN(N,2))*1 
CALL  JSCRXX<ORANCE) 

GOTO  53 

52  mXN(N. 3)  .LE.O.O)  GOTO  50 
XLOC*ALOG10(XN(N .3 ) ) ♦ 1 
CALL  JSCRXX(COLD) 

53  PO  55  NC= 1  ,4 

CALL  CSARLI (20) 

RAD=0. I2*FL0AT<NC)*XL0G^4. 

CALL  CPAR2D ( XHOPE ( N ) . YNODE (N) . RAD .0.0.6. 283 ) 

55  CONTINUE 

RAD=0.I2*XL0C 

CALL  CPAR2D <  XNODE <  N  > . YNODE ( N ) . RAD .0.0.6. 283 ) 

50  CONTINUE 

HCOUNT= 1 
COTO  I 

100  CONTINUE 

call  JaluaztO.) 

RETURN 

END 


Figure  C-3.  (Continued.) 


SUBROUTINE  DSCALC<NTH.T,TUPD.NDOSER.PrV> 

C  SUBROUTINE  TO  CALCULATE  DOSES 
C  THIS  IS  SET  UP  FOR  LLL  DATA 
C  MARCH  1982 
C 

c 

INCLUDE  ‘ DSKDiPARAM.COM' 

INCLUDE  ’DSKDiDOSE.COM’ 

INCLUDE  ’DSKD.NODE.COM' 

I  NCI.UDE  'DSKDiLlNK.COM' 

INCLUDE  ’DSKDtROAD.COM’ 

Include  ’dakd t text .con’ 

DIMENSION  P DOSE <40. 40.3> ,PPDOSE(40 .40.3) 

C  SET  UP  PAST  DOSE  MATRIX  TO  ALL  ZEROS  FIRST  TIME  IN 
1F1NTM  .NE.  DCOTO  10 

C  SET  FRAC  TO  CORRESPOND  TO  CORRECT  TIME  PERIOD 

C  PPV  IS  PEOPLE  PER  VEHICLE— ALL  XL  AND  XN  VARIABLES  ARE  IN  VEHICLES 
FRAC»T^TUPD*PPV 
DO  20  !«1 .40 
DO  20  J«I .40 
DO  20  K* 1 .3 
P DOSE (I ,J.K)-0. 

2#  PPDOSEC I , J,K)*0. 

COTO  60 


Figure  C-4.  Subroutine  DSCALC. 


10  CONTINUE 

iFlNDOSER  .NE.  1 )COTO  60 

DO  50  1*1.40 

DO  50  J*1 .40 

PDOSE! 1  .  J.  I >*D0SE< I . J . 1 ) 

DO  50  K =  2  3 

PDOSE! I . J ! K ) *  DOPE  < I .J.K'-PPDOSE! I .J.K) 

50  PPDOSE* I .J.K)*DOSE!l .J.I  > 

60  CONTI NUE 

DO  70  N* 1 , NRINCS 
RINrOFCNI'O.O 
70  CONTINUE 

C  CALCULATE  DOSE  COMPONENT  FROM  NODE  IN  AND  OUT  CUES 
C  NOTE  DOSE! I .J. I >VAS  GIVEN  AS  A  RATE  INHEMsKlt).  NOT  MR EM 
C  CORRECTION  IS  MADE  FOR  THIS 
DO  30  N* I .NNODES 
NR*N0RZ0N!N> 

NX*NGX<N) 

N Y*  NG Y  <  N  > 

XNUiT*>:N(N.2)+XN(N.3> 

XNUM2*PrV*XNUM 
DOS  I  *  PDOSE  <  NX .  NY .  1 )  *FRAC*T 
D0S2  *  PDOSE ( NX , NY . 2 ) *FRAC 
D0S3  =  PDOSE ( NX . NY . 3 ) *FRAC 
TDOS I =TDOS 1 +DOS 1 *XNUM 
TDCS2=TD0S2+D0S2*XNUM 
TDOS3=TDOS3*DOS3*XNUM 
R I RPOP ( NR ) *R I NPOP < NR ) +XNUM 
30  CONTINUE 

C  ADD  IN  DOSE  COMPONENTS  FOR  LINK  TRAVEL  AND  LINK  QUE 
DO  40  L= 1 .NLINKS 
NX=NCX(LNOUT<L) ) 

N  Y= NCY < LNOUT !  L ) ) 

NR I =  NODZON ! LN I N ! L ) ) 

NR=NODZON ( LNOUT  <  L ) ) 

POPH  *  0 . 5*XL ( L . 8 » *PPV 
DOS  I =  PDOSE ( NX . NY . 1 ) *FRAC*T 
D0S2= PDOSE ( NX , NY . 2 ) *FnAC 
D0S3 = PDOSE ( NX , NY.3 )*FRAC 
TDOS 1 =  TDOS 1 +D0S1 *XL! L . 9 ) 

TD0S2*TD0S2+D0S2*XL!  L . 9  > 

TDOS3=TDOS3* D0S3*XL! L . 9 ) 

R I RPOP C NR ) *R I NPOP < NR > +XL ( L . 9 ) *PPV«TOPH 
R I NPOP <  NR I >=RINPOP(NRI I+POPH 
C  ADD  IN  TRAVEL  COMPONENT 
RSXG*NSLD<L) 

DO  40  J  *  I , NSXG 

DFRAC* D I ST < L . J > /XLEN < L ) *FRAC 

NX-IPX(L.J) 

NY*IPY(L,J) 

DOSl * PDOSE! NX. N Y. 1 >*DFRAC*T 
D0S2=PD0SE(HX.NY ,2>*DFRAC 
D0S3*PD0SE<  NX , NY . 3  > *DFRAC 
TDOS I *TDOS I +DOS 1 *XL< L . 8 ) 

TD0S2*TD0S2+D0S2*XL<  L . B) 

TD0S3*TD0S3>  D0S3*XL (L.8) 

40  CONTINUE 

RETURN 
END 


Figure  C-4.  (Continued.) 


SUBROUTINE  ERROR 
TYPE  10 

10  FORMAT!/”  INCORRECT  RESPONSE  -  TRY  AGAIN1) 
RETURN 
END 


Figure  C-5.  Subroutine  ERROR. 


SUBROUTINE  HISTOCRINPOP. tin*) 

C  THIS  ROUTINE  PLOTS  HISTOGRAMS  OF  POPULATION  IN  ONE  MILE 
C  RINGS  AROUND  THE  SITE. 

C 

INCLUDE  * S  YS*L I BRAR Y i COLORDEF * 

INCLUDE  *DSKD: TEXT. COM' 

DIMENSION  RINPOPO  1  )  ,X<2>  .  Y<2) .  1 1 1 1  le  (25 ) 
EQUIVALENCE  (X< 1 ) ,X0> . (X(2 ) .XI ) . < Y< 1 ) . YO ) . ( Y(2> . Y1 ) 
C 

call  fsxf80(0. .1 . .0. .1 . .0. .1 . .0. . 1 . ) 
call  JnagdiCgrlt') 
e 

C  label  CRID 
C 

CALL  CSCP2D< .40.0.02) 

CALL  CPTX2D( 'DISTANCE  (MILES)*') 

CALL  JSTXUPO  .  .0.  ,0.  ) 

CALL  GSCP2D(0.065. .*) 

CALL  CPTX2D( ’POPULATION*') 
call  gscp2d(.6i,.0l) 
call  Jacrxx(gold) 
call  gptx2d( 'ewrcdl* ) 

CALL  JSTXUP<0. . 1 . .0. ) 
call  Jacrxx(whi le) 

CALL  CSCP2D( .15. .95) 

CALL  CPTX2DI 'TIME  SINCE  RELEASES') 

CALL  CSCP2D1 . 15. .92) 

ENCODE ( 15. 300.LT I TLE)  TIME 

300  FORMAT  C  4X .  ’  *  '.F4.2.'  HRS* ' ) 

CALL  CPTX2DI LTITLF.) 

coll  grscp2d  (  .  15.  .09) 

encode (18.301,ltitle)  (iftx(zzln)) 

301  format  <  '  POP  .  IN  =*'.16.'  *') 
call  gptx2d 1 1 1 1 1 le ) 

call  gacp2d < . 1 5 , .86 ) 

encode ( 1 8,302 .Itltle)  (if  1 x( zzou t ) ) 

302  format ('POP.  OUT  *’,16.’  •') 
call  gp  t  x2d  <  1 1 1 1 1  e  ) 

call  gacp2d< . IS . .83) 

f  rc 1 00* 1 00 . * f  rc 

encode ( 14 ,303 . 1 t 1 t 1 e )  frclOO 

303  forma  t ( f 5 . 1 . ‘ X  EVACD  » ' ) 
call  gp tx2d ( 1 t i t 1 e ) 
call  gscp2d( . 15, .8) 

encode <21 .304 .little)  (lf!x( tpopn) ) 

304  formal! 'TOTAL  POP.  *’,17.’  •') 
call  gp tx2d ( 1 t 1 t I e ) 

call  gacp2d( .6, .95) 

call  gp t x2d( ’ACCUMULATED  POP.*’) 

cal)  gscp2d ( . 6 . . 92 ) 

call  gp t x2d ( ’  DOSES  -  MREM* ' ) 

call  gscp2d( .6, .89) 

encode (18 .401 , 1 1 1 1 1 e )  tdoal 

401  format! *  1-131 i ’ .lpe9.2. ’  »') 

call  gptx2d ( 1 1 1 1 le ) 

call  gacp2d( .6, .86) 
encode ( 1 8,402 .Itltle)  tdoa2 

402  format ( *xe-l33t ’ . Ipe9.2 , ’  •’) 
call  gpt x2d ( 1 1 1 1 1 e  ) 

call  gscp2d( .6, .83) 
encode ( 18,403, 1 1 1 1 le)  tdoa3 


Figure  C-6.  Subroutine  HISTO. 


403  format ( '  cs- 1 37  s  * , I pe9 .2 . ’  •’) 

cal)  grp  t  *2d  ( 1  1 1  1 1  c  ) 
c 

C  gr  I  d 

c 

xmin*o.o 

XNAXM6.0 
YM I N  *  1  .0 
YMAX* I .0E+06 
c 

CALL  JSCRXX  (  GREEN ) 

CALL  CR1D1  <2. XMl.it. XMAX.YMK.YMAX) 

C 

C  PLOT  HISTOGRAMS  OIT  TO  10  MILES 
C 

CALL  JSLNSZt I00./20. ) 

YO-O.O 
DO  10  1*1 .10 

1 f trlnpopt I ) . le .0)  goto  10 
Yl =ALOC! 0<RI NPOTl I )) 

X©=  I 
XI  *X0 

1RTN  =  JPPL2A<X.Y.2> 

10  CONTINUE 

C 

C  PLOT  HISTOGRAM  AT  15  MILES  FOR  THE  REST  OF  THE  POPULATION 
C 

CALL  JSCRXX t GOLD) 

If (rlnpopt 1 1 ) . le .0.0)  goto  20 
Yl =  ALOG 1 0 ( R I NPOP  <11)1 
X0a I  5. 

XI=XO 

CALL  JPPL2A(X.Y.2) 
c 

20  CALL  JSLNSZt 0)  ! RESET  LINE  WIDTH 

C 

RETURN 

END 


Figure  C-6.  (Continued.) 


SUBROUTINE  1NITIA 

C  SUBROUTINE  TO  INITIALIZE  PLOTTING  METHODS 
C  THIS  VERSION  IS  FOR  DJSPLA  PACKACE 
C 

c 

1 NCLUDE  • SYSOL I BRARY : COLOnDEF  * 

1GRID  »  CREEN 
e  CALL  BCNPL(O) 

X0R1G=-1 3 . 

Y0R1G=-1 5 . 

XSTEP*3.75 

YSTEP*3.75 

xup  *  13. 

yup  *  13. 

xutai  *  047. 

xdtt  *  30 . 

yutai  *  3678. 

yalze  *  30. 

xtanutm  *  xu  tm 

xmxiitnj  '  xslze  ♦  xutm 

ynnuta  *  yutm 

ymxutm  «  yalze  ♦  yutnj 

call  gsxfdflO) 

call  f azf80(0. . 1 . .0. . 1 .  .0. .  1  .  ,0.  .  t  .  1 
C 

c  Initialize  grafeore 

c 

call  Jnegdof  *  ) 


Figure  C-7.  Subroutine  INITIA. 
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c  CALL  TITLE( * SAMPLE* .6. 'LAT' ,3. 'LONC' .4.0.8) 

CALL  JSCRXXi I GRID) 

call  rrldl ( 0 , zor l g , xup . yor 1 g , yup ) 

C  call  grgr2d<0.xcrtg.xup.yorig;.>ntp,:nnRntm,xmxutn.ynnutm,yRxntm! 

c  CALL  CU A PH ( XOR l G . XSTET . VOR 1 G , YSTKP 1 

c  CALL  FRAME 

return 

END 


Figure  C-7.  (Continued.) 


SUBROUTINE  ISOPLTt  tsext ) 

SUBROUTINE  TO  PLOT  APPROPRIATE  ISOCURVES 

SULROUTIRE  TAKES  CR1D  VALUES  AS  INPUT  AND  DEVELOPS  ISOCURVES 


INCLUDE  'DSKD: BOXES. CON’ 
Include  'dekdtacleve.com* 


dimension  vt<10) 

data  nx.ny.lleva  /  40.40.4/ 

READ  CR1D  INITIALIZATION  DATA 

50  CALL  PLTSETt INUC.TNEXT) 

IF< INUC.EQ.O)  CO  TO  900 

call  Jivt2d(vt) 

call  lexfOOU  .  .40.  ,1  .  .40 . .  VT<  1 )  .  v  1 1 2) .  v  t  (3 ) .  vt  <4)  > 
call  gpc tel (acleve . 1  leva ,  u,  nx.  ny) 

cal  1  faxfBOtvt  <5) ,vt(6),vt<7),vt'8>,vt(l).vt(2>.vt(3>.vt<4*> 

*>0  continue 
return 
end 


Figure  C-8.  Subroutine  ISOPLT. 


_  SUBROUTINE  LNKSET 

SUBROUTINE  TO  DETERMINE  LENGTH  OF  EACH  LINK  IN  EACH  CRID 
FEBRUARY  1982 


INCLUDE  'DSKDtPARAN.COM' 

INCLUDE  'DSKDtNODE.COM' 

INCLUDE  'DSKDtLINX.COM' 

DIMENSION  Xt<M2),X2(M2),Yl<H2).Y2(M2).XP(H2),YP<H2>,IPt<H2>, 
*  I P2 ( M2 ) 

C  LOOP  ON  LINKS 

_  DO  10  L«l .BLINKS 

C  DETERMINE  LINK  BECINNINC  AND  END  NODES.  AND  CRID  LOCATIONS 
NIN*LNIN(L> 

NOUT»LNOUT(L> 

NCXB*NCX(NIN) 

NCXE*NCX(NOUT) 

NCYB*NCY(NIN> 

NCYE*NCY(NOUT) 

XD'XNODE(NIN) 

XE«XNODE<NOUT) 

YB«YN0DE(N1N) 

YE* YNODE ( NOUT  > 


Figure  C-9.  Subroutine  LNKSET. 
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C  CALCULATE  X  AND  Y  LENGTHS  OF  LINK 
DX«XB-NE 
DY*YB-YE 

XLEN ( L ) *  SQBT ( DX*  DX+DY^DY > 

C  CALCULATE  DELTAS  FOR  GRID  LOCATIONS 
HX*NGXD-NCXE 
NY*NCYB-NCYE 

C  CALCULATE  SLOPE  AND  Y- INTERCEPT  IF  DX  IS  NOT  «  0 


IFIDX  .EQ.  0.)GCT0 
SLOPE  ^DY/'DX 
B*YB-SLOPE*XB 

19 

19 

CONTINUE 

C  SET 

UP  DIRECTIONALITY  OF 
IF! MX )  31,32.32 

SEARCH  FOR  POINTS 

31 

ND!RX*1 

MCOflX  =  0 

COTO  33 

32 

MDIRX--1 

ncop,x=  i 

33 

CONTINUE 

IEI NY )  34.35.35 

34 

PDIRY* 1 

MC0RY«0 

GOTO  36 

35 

MDIRY*-I 

MCORY- I 

30 

CONTINUE 

C  SET  UP  POINTS  INTERCEPTED  BY  X  CRID  LINES 
NP=  I 

XI (NP)=XB 
Yl (NP)=YB 
IP1 < NP )  =0 

IFl  NX  .EQ.  )GOTO  41 
NCXEM-HCXE- WDIRX 
DO  40  N  =  NGXU , NCXEN . HD I RX 
NP=NP+1 

XI <NP)*CRIDX< N-MCORX ) 

Yl ( NP ) =SLOPE*XI <  NP ) +B 
IPI I NP  >=N 

IFIXIINP)  .EO.  XB  .AND.  Y1<NP)  .EQ.  YB)NP»NP-I 

40  CONTINUE 

I  FIX  I INP)  .EQ.  XE  .AND.  YHNP)  .EQ.  YE)  GOTO  42 

41  CONTINUE 
NP»NP*T 
XIINP)«XE 
Y1INP)»YE 
IPl INP)*NGXE 

42  CONTINUE 

C  TYPE  * , (XI  IK) ,K* 1 ,NP) 

C  TYPE  ».(Y1(K).K=I.NF> 

IFIMY  .EQ.  0)G0T0  51 
MP*0 

C  SET  UP  POINTS  INTERCEPTED  BY  Y  GRID  LINES 
NCYLM  s  NG YE-N  D I RY 
DO  50  M*  NCYB . NGYEN , HDI RY 
HP*MP«-1 

Y2(NP)«CniDY( N-MCORY ) 

iP2<np>>n 

IFIDX  .EQ.  O.JCOTO  52 
X2 1 HP ) *  I Y2 1 MP ) -B ) ✓SLOPE 
COT O  50 

52  X2IMP)=XE 

50  CONTINUE 

C  TYPE  *.(X2(K).K*I .HP) 

C  TYPE  *,(Y2(K).K*l.HP) 

51  CONTINUE 

C  SORT  OUT  ALL  CRID  I NTERCF.PTS < X  AND  Y> 

NS«0 
MSTART* I 
IXELDG*0 
DO  NO  N»I,NP 
IYKLDC'0 


Figure  C-9.  (Continued.) 
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X=Xl (X) 

V* V I ( N  » 

IF<N  .EQ.  DCOTO  74 
If(MY  .EO.  0 )GOTO  71 
f1Y8rMS7ART 

1 E(NYB  .EO.  HP*DCOTO  71 
1)0  70  I1=MYB.NP.  1 
1 YKLDC -N 
XX*X2(M) 

YY=Y2  <  fl ) 

IFCXX  .LT.  XLAST  .AND.  XX  . LT .  X  OR. 

*  XX  .CT.  XLAST  .AND.  XX  .CT.  X  >GOTO  71 
NSTART=MSTART-M 

ir<XX  .EO.  XB  .AND.  YY  .EQ .  YBIOOTO  70 

NS=NSD 

XP<NS)=XX 

YP<  NS ) s  YY 

IPY<L.NS-1 ) *  I P2 ' M ) 

I  F<NS  .NE.  2 )COTO  76 
IPX<  L.NS-1  )*NGXF. 

COTO  77 

76  CONTINUE 

IFnXELDG  .EQ.  0 )COTO  86 
I PX<  L . NS-1 )  r  IP  I <  N ) 

COTO  77 

C6  CONTINUE 

IPX»L.NS-1 >=IPX<L.NS-2) 

77  CONTINUE 
XLAST=XX 
YLAST-YY 

70  CONTINUE 
I YKLDG'O 

1F<N  .NE.  NPJGOTO  71 

1F(X  .EO.  XLAST  .AND.  Y  -EO.  YLAST)COTO  60 

NS=NS+1 

XP<  NS ) =X 

YP<  NS) = Y 

IPXlL.NS-l )*NGXE 
IPY<L.NS-U»NCYE 
COTO  60 

71  CONTINUE 

IF<X  .EQ.  XLAST  .AND.  Y  .EQ.  YLAST /GOTO  60 
1F<  1 YKLDC  .NE.  OIXKLDCM 
74  CONTINUE 

NS=NS*I 
XP(NS) =X 
YP(NS)=Y 

C  1F(N  .NE.  DCOTO  72 

C  I PX  <  L , NS ) =NGXB 

C  IPY<L.NS)*NCYB 

C  COTO  73 

72  CONTINUE 

IFINS  .EQ.  1  .AND.  N  .NE.  NPICOTO  73 
I PX< L.NS-1 )*IP1 (N) 

I  FINS  .NE.  2  >GOTO  78 
I F *  1 YKLDC  .NE.  8>C0T0  95 
IPY<  L.NS-1 )*NGYB 
COTO  73 

78  CONTINUE 

1 F  <  I YKLDC  .EQ.  0)C0T0  79 
♦5  I PY« L.NS-1 )«IP2< I YKLDC) 

C  IXKLDCM 

COTO  83 

79  CONTINUE 

IFCMSTART  .NE.  HP*- DCOTO  94 
I PY I L.NS-1 )*NCYE 
COTO  96 

94  CONTINUE 

I PY<  L.NS-1 >* IPY1L.NS-2) 

96  CONTINUE 

IXKLDC'O 
83  CONTINUE 


Figure  C-9.  (Continued.) 
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CONTINUE 
XLAST'X 
YLAST* Y 
6©  CONTINUE 

1 F (MX  .EQ.  0  .AND.  MY  .EQ.  OHPX'L. 1 >=NCXB 
C  DETERMINE  NUMBER  OF  SEGMENTS  PER  LI NK! NSLD)  . LENGTH  OF  EACH 
C  SEGMENT) DI ST) ,  AND  GRID  CELL  LOCATION ( NGLOCX  AND  NCLOCY) 

C  TYPE  *,<XP<JO  .K* 1 .NS) 

C  TYPE  *. (YP(K> .K' 1 ,NS) 

NC*0 

DO  80  N*  1 .NS 
IFlN  .EQ.  1 )COTO  81 
NCSNC*T 
XA=XP<  N ) 

YA*YP<N> 

XDIF*XA-XZ 

YDIE=YA-YZ 

DIST!  L .  NC ) =SQRT  <  XDI F*XDIF+YP J  F*YDI F ) 

I F( PI  ST ( L . NC )  .NE.  0.)C0T0  82 
NC=NC- I 

82  CONTINUE 

81  CONTINUE 

XZ=XP<  N  > 

YZ*  YP<  N  > 

80  CONTINUE 

NSLDIL)*NC 
LTOT  *  NC 

WRITE! 1 3 , 2000 ) L . NSLD ( L ) 

DO  90  N=  1  .LTOT 

*0  WRITE! 1 3 , 260 1)1 PX ( L . N ) . IPY(L.N) .DISTIL. N) 

2001  F0R(1AT(2I8.  I  PE  1 0.3) 

2000  FORMAT! 16, 16 ,/lCI 6. '5E10 .4) 

10  CONTINUE 

CLOSE!UNlT=l3) 

RETURN 

END 


Figure  C-9.  (Continued.) 


SUBROUTINE  RODSET 

C  SUBROUTINE  TO  DETERMINE  WHERE  NODES  ARE  LOCATED  RELATIVE  TO  CRID 
C  FEBRUARY  1982 
C 
C 

INCLUDE  ’ DSKD: PARAM.COM ’ 

INCLUDE  'DSKD: NODE. COM’ 

OPEN! UN IT* 13 . f I le=  *4ak  tEVACDB.NLD* ,.t»tn»* 'new' ) 

C  SET  UP  CRID  SPACING 

DO  40  N= I , NCRIDX 

40  GR I DX  <  N ) 1 DELX* FLOAT ! N ) +XCORN 

DO  50  H* I , NCR  1 DY 

50  CRI  DY  ( N )  =DF.LY*FLOAT(N )  +YCORN 

C  LOOP  ON  NODES 

DO  10  N* 1 .NNODES 

C  CALCULATE  X  AND  Y  DISTANCE  FROM  ZERO  REFERENCE  FOR  EVACUATION 
DX« XNODE < N ) -XZERO . 

DY* YNODE I N ) -YZERO 
RAD* SORT 1 DX*DX+PY*DY ) 

C  DETERMINE  WHICH  SECTOR! 1-16)  NODE  IS  LOCATED  IN 
Z* ATAN2 ! DX , DY ) * 57 . 296 
IF!Z  .LT.  0. >Z«Z+360. 

Z*Z*I I .25 

IFIZ  .CE.  360. )Z*Z-360. 

NODSEC ! N ) *  I F I X ( 22 .  5 )  *  1 
C  DETERMINE  WHICH  RING! 1-5)  NODE  IS  LOCATED  IN 
DO  20  K* I .NRINCS 
IF! RAD  .CT.  RINC!K))COTO  20 
NODZON!N)*K 
GOTO  21 


Figure  C-10.  Subroutine  NODSET. 
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20  CONTI  NUE 

Typr.  looo.w 

1000  FORMAT! *  RODE  RUHBER  ’.13.’  IS  OUT  OF  RANGE’  > 

NERKORM 

21  _  CONTINUE 

C  DETERMINE  WHICH  NODES  ARE  TO  BE  EVACUATED 
C  CONSIDER  ONLY  A-TYPE  NODES 

1F(N0D1D(N)  .NE.  ‘A  • ) GOTO  32 
DO  30  J- 1 , NZONES 
NS=NSEC< J ) 

IF<RAD  .CT.  ZONE ( J ) )COTO  30 
DO  3 1  I  =  I . NS 

I F  ( NODSEC  ( N )  .NE.  MSECS (J . I))COTO  31 

NODVAC  <  N  )  s  1 
COTO  32 

31  CONTINUE 

30  CONTINUE 

32  CONTINUE 

C  SET  UP  SQUARE  GRID  INDICES 
DX=XNODE<  N ) -XC0KN 
DY  *  VNODE ( N ) - YCOPN 
NCX( N ) -  I  FI X( DX'DELX) +1 
NCY <  N  )  =  I  F I  X<  DY/DELY  >♦ 1 

WRITE! 1 3 . 2000 ) N . NOD I D ( N ‘ , NODSEC ( N ) . NODZON ( N ) .NODVACCN) . 
*  NCX(N) .NCY(N) 

2000  FORMAT! 16. A! .516) 

10  CONTINUE 

RETURN 
END 


Figure  C-10.  (Continued.) 


SUBROUTINE  PLTEND 

CALL  JESCDAl 'GRID' ) 

CALL  JXCLFR 

CALL  JESCXXl * CR 1 D ' ) 

RETURN 

END 


Figure  C-ll.  Subroutine  PLTEND. 


SUBROUTINE  PLTSETUNUC  .TNEXT) 

C  SUBROUTINE  TO  SELECT  AND  PLOT  CURVES 
C 

c 

INCLUDE  ’SYS* LIBRARY : COLORDEF ’ 

INCLUDE  ’DSKDtBOXES.COM’ 

INCLUDE  ’DSKDtDOSE.COM’ 

Include  *dskd : eel  eve . con’ 

Include  *dskd t text .con’ 
dimension  aclev<4,3)  !•» 

DIMENSION  LT1TLE120) 

C 

read  (8.*)  nclev  tread  contour  levels 


c 

C 

C  FIND  OUT  WHAT  TYPE  OF  ISOPLETH  SHOULD  DE  PLOTTED 
C 


10  TYPE  1001 
lOOl  FORMAT </V*  ENTER 
2/"  ENTER  •I"  FOR 
3/ ’  ENTER  FOR 
4/’  ENTER  "S"  FOR 
5/’  ENTER  "©•  FOR 


NUCLIDE  TO  BE  PLOTTED. ’✓ 
1-131 .  OR’ 

XE-133.  OR’ 

CS- 137 ,  OH’ 

NO  PLOT  AT  THIS  TIME  STEP.’) 


8* 


Figure  C-12.  Subroutine  PLTSET. 
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o««  m  •*  oftft  non 


ACCEPT  s.INUC 
I F 1 1 NUC .  EQ  .  0 )  RETURN 

I F 1 1 NUC . EQ . 1 .  OR  .  I NUC  .  EQ  .  2 .  OR .  I  NUC  .  EQ .  3 )  GOTO  13 
CALL  ERROR 
CO  TO  10 

LOAD  THE  APPROPRIATE  VALUES  FOR  THIS  PLOT  INTO  ARRAY  U. 

IS  DO  20  J  =  l .40 
DO  20  1=1. 40 

20  U! 1 ,J)*DOSE! 1 . J. INUC) 
do  16  i=l,l levs  ISO 

aclevell)  =  sc  1 ev ( 1  , 1 nuc )  !** 

1 6  cont 1 nue  ! *1 

SET  UP  THE  TEXT 

CALL  JSCRXX (GREEN) 

CALL  CSCP2D<6. . 14. ) 

IF< INUC.EQ. 1 )  CALL  CPTX2D( *1-131  INHALATIONS*) 

IF< 1NUC.EQ.2)  CALL  GPTX2D( 'XE-133  DIRECTS') 

IF! 1NUC.EQ.3)  CALL  GPTX2D! *CS- 1 37  INHALATIONS') 
CALL  CSCP2D< 6 . . 13. ) 

CALL  CPTX2D( ’DOSE  ISOPLETHSS*) 

CALL  CSCP2D(6. . 12. ) 

F.NCODE<  12,150.  LTITLE)  TNEXT 
50  FORMAT! 'AT' .F5.2. *  HRSS') 

CALL  CPTX2D  ( LT 1 TLF. ) 

CALL  JSCRXX ( MAGENTA ) 

CALL  CSCP2D! 10. . 1 1 . ) 

ENCODE! 9 .200 . LTITLE)  ACLEVS! 1 ) 

CALL  CPTX2D! LTITLE) 

CALL  JSCRXX !LT  MAGENTA ) 

CALL  GSCP2D! 10. . 10. ) 

ENCODE (9. 200. LTITLE)  ACLEVS (2) 

CALL  CPTX2D! LTITLE) 

CALL  JSCRXX ! MED_MACENTA ) 

CALL  C-SCP2D!  10.  .9.  ) 

ENCODE 19.200. LTITLE)  ACLEVS !3) 

CALL  GPTX2D! LTITLE) 

CALL  JSCRXX!  MI  I TE) 

CALL  CSCP2D! 10. .0. ) 

ENCODE 19. 200, LTITLE)  ACLEVS 1 4 ) 

CALL  CPTX2D! LTITLE) 

SO  FORMAT! 1PE8. 2. *S* ) 

RETURN 

END 


Figure  C-12.  (Continued.) 


SUBROUTINE  SITE! ISITRD, TIME) 
SUBROUTINE  TO  PLOT  OUT  SITE  SPECIFICS 


INCLUDE  *8YSSLI BRARY i COLORDEF ’ 

INCLUDE  • DSKD.PARAM.COM* 

INCLUDE  ’DSXD: NODE. COM* 

INCLUDE  ’DSHD : LINX.COM* 

Include  ’dskditezt.com* 

DIMENSION  XSITE16.46) ,YSITE(«,46) .NSITE16) , 1 tltle<23> 
DIMENSION  XPLT146) ,YPLT<46) .xxx<2) .yyy(2) 

equivalence  (xxx( 1 ) ,xl ) , !yyy( 1 ) ,yl ) , ! xxx(2) ,x2), (yyy<2) ,y2) 

DATA  NSITE/13, 41 ,19.3,3. I5x 

I ROAD  *  YELLOW 

I  LAKE* CYAN 

ICONTY«VHITE 

I SITE* ORANGE 

ICRI D*CREEN 


Figure  C-13.  Subroutine  SITE. 
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Cl  NR* 10.0 
COUTR* 15.0 

OPEN (UN  IT* 18, f 1 le* 'DSKrSITE. INP  ’ .status* ’old' ) 
IF ( 1SITRD.NE.0)  CO  TO  30 
C  READ  SITE  DATA 
«  CALL  COMPRS 

« 

call  gndlmdlO) 

«  call  vailld  'COLOR  MONITOR 

CALL  LV11ID  t HARDCOPY 
call  (*rtU<0)  tno  label*  on  contour* 

DO  10  1*1,6 

DO  10  J* 1 ,NSITE( I ) 

10  READ<  1Q.1000>  K.XSITEU  .J)  .YSITEd  .J) 

CLOSE (UN  IT* 18) 

C  PLOT  THE  SITE  DATA 

CALL  JSCRXX! ICR1 D) 

30  CALL  INJTIA 

call  £sarll(200>  (circle  segment  count 

call  gpnr2d (0.0.0.0.CINR.0.0.6. 203 > 
call  gpar2d (0.0. 0.0. COUTR. 0.0. 6. 200) 

CALL  JSCRXX( ISITE) 

call  gsarli(20)  (circle  segment  count 
call  gpar2d (0 .0,0. 0.0. 3. 0.0. 6. 283 ) 

DO  40  1*1.6 

IFO.EQ.2)  CALL  JSCRXX  (  I  LAKE ) 

IFO.GT.2)  CALL  JSChXX <  1 CONTY > 

DO  50  J*1 .NSITE( I ) 

XPLT( J ) =XS1TE< 1 .J) 

50  YPLT(J)*YSITE( I . J) 

call  gpcv2d ( xpl t ,ypl t ,ns l te ( i ) ) 

♦0  coni  1 nue 

C  PLOT  OUT  HOAD  NETWORK 

DO  100  L* 1 .RLINKS 
N I N  =  LN I N  ( L ) 

NOUT*  LNOUT ( L  > 

XI*XNODE(NIN) 

YI=YNODE(NIN) 

X2=XN0DE(N01T) 

Y2* YNODE ( NOUT) 

c  CALL  RLVEC( XI .Y1 ,X2.Y2, 0000) 

CALL  JSCRXX  < I  ROAD ) 

CALL  CSARLI ( B> 

CALL  CPAR2D(X1 ,Y1 ,0.12.0.0.6.233) 

Irtn  =  Jpp 1 2a ( xxx , yyy . 2 • 

CALL  CPAR2D(X2,Y2.0. 12.0.0,6.263) 

100  CONTINUE 

1000  FORMAT (14. 2F 10.4) 

1 f ( 1st  t rd .gt .0)  CO  TO  90 
call  Jacrxx( Igrid) 
call  gscp2d(-14. , 14. ) 
encode (24. 200. 1 1 1  tie )  tno tic 

200  format ( 'WARN INC  TIME  *',r5.2.'  HRS*  * ) 
call  gptx2d ( 1 t 1 t 1 e ) 

cal)  gscp2d ( -14 . , 13 . ) 
encode (21 ,201 . 1 t I t le )  tdelay 

201  format( ’PREP  TIME  *',F5.2,’  HRS* ’ > 
call  gptx2d< 1 ti t le) 

call  gacp2d (-14. ,12.) 

encode (21 , 202 .ltltle)  (iflx(tpopn)) 

202  format (’TOTAL  POP.  *’,17.'  •’) 
call  gptx2d( 1 tl tie) 

RETURN 

*0  CALL  JSCRXX (CREEN) 

CALL  CSCP2D( -14 . ,14.) 

CALL  CPTX2D1 'TIME  SINCE  RELEASE*’) 

CALL  CSCP2D(-14. .13. ) 

ENCODE (15, 300  , LT ITLF. )  TIME 
300  FORMAT <4X , ’ *  ’.F4.2.’  HRS* ' ) 

CALL  CPTX2D( LTITLE) 

call  gscp2d ( -14. , 12 . ) 

encode ( 17 ,301 .ltltle)  (iflx(zzln)) 


Figure  C-13.  (Continued.) 
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301  formate  ’POP.  IN  ■*.!«, *•*) 

cal)  grptx2d<l  11  tie) 

call  *scp2<l<-14.  .  1 1 .  ) 

*ncod» ( 17 ,302, Itltle)  (I f I x( zzou  t ) ) 

302  formate ‘POP.  OUT  ■,.I6.*0,> 

call  gptx2d( 1 t l t le) 

call  *»cp2d(-I4. , 10. ) 

fre I00=l00-*frc 

encode ( 14 ,303. 1 1 1  * le )  fre 100 

303  format efS. 1 . EVACD  •’) 
cal)  gptx2d ( 1 1 1 1 le ) 

RETURN 

END 


Figure  C-13.  (Continued.) 


C  ACLEVS.COM 

COMMON  'ACLEVS  ACLEVS < 4 > . I LEVS 


Figure  C-14.  Common  tile  ACLEVS.COM. 


C  BOXES.COM 

COMMON  XBOXES'  XCRIDC40) ,YCRID(40) , XORIC , YORIC .XDELT 
t ,YDELT,U<40,49) 


Figure  C-15.  Common  file  BOXES.COM. 


C  DOSE.CON 

COHMON'DOSEXDOSE C 40 . 40 . 3 ) 


Figure  C-16.  Common  file  DOSE.COM. 

| 

|  C  LINK.COM 

!■  C  COMMON  FOR  LINK  DATA 

COMMON'LNKl/NLlNKS,LNINCN2)  ,LN0UT(N2)  ,DIST(N2,H1 ) ,  IPX(N2.M1 )  , 
.  «  IPY<N2,M1 ) ,NSLD(N2) ,LNUML(N2) ,xlen(n2) 


Figure  C-17.  Common  file  LINK.COM. 


C  NODE.COM 

C  COMMON  FOR  NODE  DATA 

COMHON'NODl XNNODES , XZEHO , YZERO , XNODE (N 1 ) , YNODEeNl ) ,N0DSEC(N1 >  , 

*  RINCtNS) , NODZON ( N 1 > .NODIDCNl > .NSECe NS) .ZONEINS) . 

*  MSECS ( NS, N4 ) ,N0DVAC(N1 ) .NHINCS.NZONES.NERROR. 

*  NCR I DX , NCR  1 DY , XCORN , YCORN . DELX , DELY . NCX1N 1 ) , NCYCNl > , 

*  CR1DX(N6> ,CRIDY(N6) .rlnpop(nS) 


Figure  C-18.  Common  file  NODE.COM. 
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PAHAHETDI  Kl«209,R2»2©d,H3*29,R4*16,H5M3,H6»4#,Hl»2B,R2»W 


Figure  C-19.  Parameter. 


C  ROAD. COM 

G0HMHI/1MNUI/XL(1I2, 12)  ,XR(llt  .3) 

Figure  C-20.  Common  file  ROAD.COM. 


e  text. com 

c  cobbou  for  text  for  dleplays 
O  aleo  contain*  doaea 

couob  Xtxtl/  zzln,zzout ,  frc  ,  taot  lc ,  tdelay ,  tpopn,  tdoal 
*  tdoa2.tdoa3 


Figure  C-21.  Common  file  TEXT.COM. 
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Appendix  D: 

Input  Files  of  Sample  Case 


800. 

30. 

3. 

.3 

3. 

.3  . 

300. 

39. 

3. 

.3 

10. 

1. 

300. 

30. 

3. 

.3 

10 

.  1. 

300. 

30. 

3. 

.3 

10 

.  1. 

300. 

30. 

3. 

.3 

10 

.  1. 

300. 

30. 

3. 

.3 

30. 

3. 

300. 

30. 

3. 

.3 

30. 

3. 

300. 

30. 

3. 

.3 

30. 

3. 

309. 

30. 

3. 

.3 

30. 

3. 

300. 

30. 

3. 

.3 

30. 

3. 

100. 

10. 

1. 

.1 

30. 

3. 

100. 

10. 

1 . 

.1 

39. 

3. 

30. 

3.  . 

3  . 

03  30. 

3.  . 

30. 

3.  . 

3  . 

03  30. 

3.  . 

to. 

1.  . 

1  . 

01  30. 

3.  . 

10. 

1.  . 

1  . 

01 

30. 

3.  . 

3.  . 

3  .03  . 

003 

30 

.  3. 

.3  . 

03  . 

003 

.0003 

30. 

03  .903  .003  .0003  .00003  .000003 
.1  .91  .01  .001  .0001  .00001 
.1  .01  .01  .001  .0001  .00001 
.1  .01  .01  .001  .0001  .00001 
.1  .01  .03  .003  .0003  .00003 
.3  .03  .93  .093  .0003  .00003 
.3  .03  .93  .003  .0003  .00003 
.3  .03  .03  .003  .0003  .00003 
.3  .03  .03  .003  .0003  .00003 
.3  .03  .03  .003  .0003  .00003 
.3  .03  .03  .003  .0003  .00003 
.3  .03  .03  .003  .0003  .00003 
3  .03  .03  .003  .0003  .00003 
3  .93  .03  .003  .0003  .00003 
3  .03  .03  .003  .0003  .00003 
3  .03  .03  .003  .0003  .00003 
.3  .03  .03  .003  .0003  .00003 
3.  .3  .03  .03  .003  .0003  .00003 


Figure  D-l.  Common  input  file  ACLEVS.DAT -contour  levels  to  be  plotted  each  time  step. 


BAHCHO  SECO  TEST  CASE  3 — BUCKLEY — 01MARCH82 
3  .004  .05  .25  .85  1 

100  103  123  30 


IB 

-2. 

16 

-0. 

57 

2B 

-2. 

13 

-1. 

68 

3A 

-2. 

82 

-1 . 

71 

4A 

-2. 

83 

-2. 

75 

5B 

-2. 

84 

-3. 

86 

6A 

-2. 

85 

-5. 

99 

7B 

-2. 

73 

-6. 

88 

BA 

-2. 

70 

-a. 

12 

9B 

-2. 

70 

-9. 

24 

10Z 

-2. 

66- 

10. 

28 

UB 

-2. 

67 

-0. 

34 

12A 

-2. 

66 

1. 

15 

13A 

-2. 

69 

3. 

57 

14B 

-3. 

.23 

5. 

,17 

15B 

-4. 

.46 

6. 

,96 

16B 

-5. 

,95 

-3. 

,89 

17A 

-5. 

.96 

-4. 

.99 

18B 

-5, 

.99 

-6. 

.04 

19B 

-5. 

.97 

-6, 

.56 

20B 

-5. 

.99 

-a. 

.15 

21Z 

-5, 

.97 

-9, 

.32 

22B 

-5 

.95 

-2 

.87 

23B 

-5 

.96 

-1 

.78 

24B 

-5 

.86 

-0 

.11 

25A 

-5 

.86 

2 

.60 

26B 

-5 

.92 

4 

.70 

27B 

-7 

.05 

4 

.79 

28A 

-0 

.49 

-9 

.23 

29B 

0 

.89 

-9 

.23 

30B 

1 

.93 

-9 

.74 

3IB 

-3 

.81 

-9 

.26 

32B 

-4 

.88 

-9 

.27 

33A 

-0 

.51 

-8 

.13 

34B 

1 

.91 

-8 

.10 

35B 

-3 

.79 

-8 

.14 

36B 

-4 

.88 

-8 

.13 

Figure  D-2.  Common  input  file  EVACDB.INP  -  road  network  information. 
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37Z 

-9. 

04 

-8. 

16 

38Z 

-9. 

87 

-6. 

.72 

39B 

-7. 

66 

-6. 

.08 

40B 

-9. 

24 

-6. 

.10 

41Z-10. 

.16 

-6. 

.11 

42A 

-7. 

.64 

-4. 

,99 

43A 

-3, 

.93 

-0. 

.07 

44A 

-7. 

.63 

-0, 

.15 

43B 

-7. 

.59 

-0, 

.65 

46B-10. 

15 

-0. 

.61 

47B-10. 

.14 

-1 , 

,47 

48A 

1  . 

.10 

-3, 

,76 

49A 

-0. 

.02 

-3, 

.78 

BOA 

-1 . 

.56 

-3, 

.82 

51A 

-7. 

.42 

-3 

.91 

32B 

-9, 

.25 

-3 

.96 

53Z- 

■11. 

.00 

-3 

.97 

54A 

2 , 

.95 

-2 

.88 

95A 

0 

.56 

-0 

.51 

56B 

-6 

.93 

-1 

.42 

87Z- 1 1 . 

.84 

-1 

.46 

S8A 

-3 

.14 

6 

.05 

99B 

-3 

.93 

7 

.36 

60A 

-3 

.47 

8 

.52 

6 1 B 

-2 

.93 

9 

.73 

62B 

-3 

.73 

5 

.72 

63  B 

-6 

.30 

5 

.35 

*4  6.39O9E+02  t.eoooE+oa  1 .ooooe+oo 

67  1 .90O0E+64  l .0000E+03  I . OOOOE+OO 

98  6 . 3999E+02  1 .0O00E+03  1 .OOOOE+OO 

60  2 . 20O0E+02  t . 0000E+03  1 . 0O0OE+00 


1 

1.10 

1 

49.0 

2000. 

1 .0 

1.9 

1 .0  -1  .0 

1 

2 

0. 

0. 

2 

0.70 

1 

45.0 

2000. 

1.0 

1 .0 

1 .0  -1 .0 

2 

3 

0. 

0. 

3 

1.00 

1 

45.0 

2009. 

1 .0 

1 .9 

1.0  -1.0 

3 

4 

9. 

0. 

4 

1.10 

1 

45.0 

59. 

1.0 

1  .0 

1 .0  -1  .0 

4 

3 

0. 

0. 

9 

2.10 

1 

45.0 

2009. 

1 .0 

1.0 

1  .0  -1 .0 

5 

6 

0. 

0. 

6 

0.90 

1 

49.0 

2000. 

1  .0 

1.0 

1  .0  -1  .0 

6 

7 

0. 

0. 

7 

1.20 

1 

49.0 

50. 

1.0 

1 .0 

1 .0  -1 .0 

7 

a 

0. 

0. 

8 

1.10 

1 

43.0 

2009. 

1.0 

1.0 

1 .0  -1 .0 

a 

9 

0. 

0. 

9 

1.00 

1 

43.0 

2000. 

1.0 

1 .0 

1.0  -1.0 

9 

10 

0. 

0. 

10 

0.60 

1 

49.0 

2000. 

1.0 

1.0 

1.0  -1.0 

i 

11 

0. 

0. 

11 

1.50 

1 

45.0 

2000. 

1 .0 

1.9 

1  .0  -1  .0 

u 

12 

0. 

0. 

12 

2.30 

1 

43.0 

2000. 

1.0 

1 .0 

1.0  -1.0 

12 

13 

0. 

0. 

13 

1.70 

1 

45.0 

2000. 

1  .0 

1.9 

1 .0  -1 .0 

13 
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Figure  D-3.  Common  input  file  EVACDB.CAS  —  specific  evacuation  case  information. 


Figure  D-4.  Common  input  file  I131.LLNL-  40  x  40  doses  due  to  mI  each  time  step. 
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Figure  D-6.  Common  input  file  CS137.LLNL-40  X  40  doses  due  to  ,37Cs  each  time  step, 
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Appendix  E: 

Discussion  of  Sample  Case 

Appendix  E  discusses  the  sample  problem  presented  in  this  report.  Figure  E-l  illustrates  the  road 
network  surrounding  the  Rancho  Seco  site.  The  numbers  represent  the  103  nodes.  Table  E-l  lists  all  the 
nodes  whe>-e  people  enter  the  evacuation  road  network,  the  characteristics  of  each  node  including  x-  and 


Figure  E-l.  Initial  road  network,  with  a  node-numbering  scheme,  surrounding  the  Rancho  Seco 


y-location,  the  population  entering  the  node,  and  the  node  capacity  for  entering  vehicles.  The  total  popu¬ 
lation  involved  is  10,304.  Each  node  is  identified  as  being  an  input,  through,  or  output  node.  Input  nodes 
are  nodes  that  the  evacuating  population  enters,  and  output  nodes  are  nodes  to  which  the  evacuating 
population  travels.  The  through  nodes  are  those  through  which  the  evacuating  population  moves  during 
travel  from  an  input  to  an  output  node. 

Table  E-2  illustrates  the  characteristics  of  each  link,  including  length,  free-flow  speed,  capacity,  num¬ 
ber  of  lanes,  and  beginning  and  ending  nodes.  Potential  road  blocks  are  also  indicated.  Links  7,  55,  and  85 
have  low  road  capacities.  Since  capacity  is  the  number  of  vehicles  that  can  pass  through  a  link  per  hour, 
low  capacities  will  cause  roadblocks  during  heavy  traffic  periods.  The  low  capacities  of  links  7,  55,  and  85 
during  the  evacuation  result  in  queues  that  form  on  those  links  or  on  the  feeder  links  leading  to  them. 

The  specific  parameters  for  the  sample  problem  are  given  in  Table  E-3.  The  problem  was  run  using 
0.01-h  time  steps  and  the  results  were  reported  every  0.25  h.  The  evacuation  zone  was  up  to  10  mi  and 
included  all  sixteen  22.5°  sectors  (1  to  16). 


Table  E-l.  Node  characteristics  of  the  sample  problem  for  the  Rancho  Seco  Nuclear  Power  Plant. 

LOCATION  POPULATION 

(MILES)  ENTERING  NODE  CAPACITY 

NODE  ID  TYPE  X  Y  NODE (PEOPLE)  (VEHICLES/HR) 

g*************************** ****** ******** ***** ********************* 


1 

B 

-2.16 

-0.57 

0. 

0. 

2 

B 

-2. 13 

-1 .68 

0. 

0. 

3 

A 

-2.82 

-1  .71 

142. 

1000  . 

4 

A 

-2.03 

-2.75 

1000. 

1000. 

5 

B 

-2.04 

-3.86 

0. 

0. 

6 

A 

-2.05 

-5.99 

300. 

1000. 

7 

B 

-2.73 

-6.80 

0. 

0. 

0 

A 

-2.70 

-0.12 

409. 

1000. 

9 

B 

-2.io 

-9.24 

0. 

0. 

10 

Z 

-2.66 

-10.20 

0. 

0. 

11 

B 

-2.67 

-0.34 

0. 

0. 

12 

A 

-2.66 

1.15 

65. 

1000. 

13 

A 

-2.69 

3.57 

90. 

1000. 

14 

B 

-3.23 

5.17 

0. 

0. 

IS 

B 

-4.46 

6.06 

0. 

0. 

16 

B 

-5.95 

-3.89 

0. 

0. 

17 

A 

-5.96 

-4.99 

1000. 

1000. 

10 

B 

-5.99 

-6.04 

0. 

0. 

19 

B 

-5.97 

-6.56 

0. 

0. 

20 

B 

-5.99 

-0.15 

0. 

0. 

21 

Z 

-5.97 

-9.32 

0. 

0. 

22 

B 

-5.95 

-2.87 

0. 

0. 

-23 

B 

-5.96 

-1  .78 

0. 

0. 

24 

B 

-5.06 

-0.11 

0.. 

0. 

23 

A 

-5.06 

2.60 

600. 

1000. 

26 

B 

-5.92 

4.70 

0. 

0. 

27 

B 

-7.05 

4.79 

0. 

0. 

20 

A 

-0.49 

-9.23 

290. 

leoo. 

29 

B 

0.09 

-9 . 23 

0. 

0. 

30 

B 

1  .93 

-9.74 

0. 

0. 

31 

B 

-3.01 

-9.26 

0. 

0. 

32 

B 

-4.00 

-9.27 

0. 

0. 

33 

A 

-0.51 

-0.13 

100. 

1000. 

31 

0 

1.91 

-0.  10 

0. 

0. 

30 

B 

-3.79 

-0.14 

0. 

0. 

30 

B 

-4.00 

-8.13 

0. 

0. 

37 

Z 

-9.04 

-8.16 

0. 

0. 

30 

Z 

-9.07 

-6.72 

0. 

0. 

3V 

B 

-7.66 

-6.00 

0. 

0. 

40 

0 

-9.24 

-6.10 

0. 

0. 

41 

z 

-10.16 

-6.11 

0. 

0. 

42 

A 

-7.64 

-4.99 

345. 

tooo. 

43 

A 

-3.93 

-0.07 

152. 

1000. 

44 

A 

-7 . 63 

-0.15 

236. 

1000. 

45 

B 

-7.09 

-0.65 

0. 

0. 

46 

B 

-10. in 

-0.61 

0. 

0. 

47 

B 

-10.14 

-1  .47 

0. 

0. 
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Table  E-l.  (Continued.) 


LOCATION  POPULATION 

(MILES)  ENTER INC  RODE  CAPACITY 

RODE  ID  TYPE  X  Y  RODE  (PEOPLE)  (VEHICLES /HR) 

******************************************************************** 


40 

A 

1  .  to 

-3.76 

22. 

500. 

49 

A 

-0.02 

-3.70 

16. 

500. 

50 

A 

-1  .56 

-3.02 

1000. 

500. 

51 

A 

-7.42 

-3.91 

345. 

1000. 

52 

B 

-9.23 

-3.96 

0. 

0. 

53 

Z 

-1 1 .00 

-3.97 

0. 

0. 

54 

A 

2.93 

-2.00 

1 1  . 

500. 

55 

A 

0.56 

-0.51 

1000. 

1000. 

56 

B 

-6.93 

-1  .42 

0. 

0. 

57 

Z 

-11.04 

-1  .46 

0. 

0. 

50 

A 

-5.14 

6.05 

639. 

1000. 

54 

B 

-3.93 

7.36 

0. 

0. 

6ft 

A 

-3.47 

0.52 

220. 

1000. 

6! 

B 

-2.93 

9.73 

0. 

0. 

62 

B 

-5.73- 

5.72 

0. 

0. 

63 

B 

-6.30 

5.35 

0. 

0. 

64 

A 

-0.62 

3.63 

639. 

1000. 

65 

B 

-9.83 

2.63 

0. 

0. 

66 

Z 

-12.43 

-0.13 

0. 

0. 

67 

A 

-5.01 

4.77 

1000. 

1000. 

60 

B 

-4.83 

5.18 

0. 

0. 

69 

B 

-7.54 

5.42 

0. 

0. 

70 

B 

-0.61 

4.91 

0. 

0. 

71 

Z 

-10.39 

6.36 

0. 

0. 

72 

A 

0.22 

9.01 

12. 

1000. 

73 

B 

-2.07 

9.17 

0. 

0. 

74 

Z 

-4.27 

11.19 

0. 

0. 

75 

B 

2.47 

7.49 

0. 

0. 

76 

B 

3.82 

6.33 

0. 

0. 

77 

A 

4.89 

5.41 

8. 

1000. 

70 

B 

7.84 

4.06 

0. 

0. 

79 

B 

8.05 

2.86 

0. 

0. 

OO 

B 

-0.96 

0.52 

0. 

0. 

01 

B 

0.44 

0.53 

0. 

0. 

02 

B 

2.03 

1.16 

0. 

0. 

03 

A 

3.96 

3.15 

0. 

1000. 

64 

A 

7.59 

3.32 

50. 

1000. 

85 

B 

11.30 

1  .69 

0. 

0. 

CO 

Z 

12.03 

0.77 

0. 

0. 

87 

B 

-3.10 

-0.84 

0. 

0. 

80 

B 

4.30 

0.50 

0. 

0. 

09 

Z 

4.45 

10.22 

0. 

0. 

90 

A 

9.13 

-0.89 

500. 

1000. 

91 

B 

10.10 

-0.13 

0. 

0. 

92 

B 

10.49 

1  .03 

0. 

0. 

93 

A 

9.14 

-2.03 

30. 

1000. 

94 

B 

9.75 

-2.41 

0. 

0. 

95 

B 

11.73 

-1.55 

0. 

0. 

96 

B 

6.93 

-4.47 

0. 

0. 

97 

B 

4.57 

-6.77 

0. 

0. 

90 

B 

2.99 

-9.80 

0. 

0. 

99 

A 

1  .92 

-7.29 

75. 

1000. 

100 

Z 

1  .66 

-10.31 

0. 

0. 

101 

B 

-0.47 

-10.25 

0. 

0. 

102 

B 

0.92 

-9.19 

0. 

0. 

103 

B 

0.94 

-10.30 

0. 

0. 

RODE  TYPE i  A  *  INPUT  NODE  RELEASE  POINT  IS  AT  X»0  AND  Y*0 

B  «  THRU  RODE 
Z  *  OUTPUT  NODE 
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Table  E-2.  Link  characteristics  of  the  sample  problem  for  the  Rancho  Seco  Nuclear  Power  Plant. 


LENGTH  FREE  FLOW  CAPACITY  START  END  NUMBER  OF 

LINK  ID  <M1LE)  SPEED (MPH )  (VEIIICLES^HR)  RODE  NODE  LANES 

**********il^:|^***«*******x****K***************»****«*********s***********x 


t 

1.10 

45. 

2000. 

1 

2 

2 

0.70 

45. 

2000. 

2 

3 

3 

1 .00 

45. 

2000. 

3 

4 

4 

1.10 

45. 

1000. 

4 

5 

5 

2.10 

45. 

2000. 

3 

6 

6 

0.90 

45. 

2000. 

6 

7 

7 

1  .20 

45. 

150. 

7 

8 

8 

1.10 

45. 

2000. 

8 

9 

0 

1  .00 

45. 

2000. 

9 

10 

10 

0.60 

45. 

2000. 

1 

11 

n 

1.50 

45. 

2000. 

11 

12 

12 

2.30 

45. 

2000. 

12 

13 

13 

1  .70 

45. 

2000. 

13 

14 

14 

1 .50 

45. 

2000. 

14 

15 

15 

1  . 10 

45. 

2000. 

16 

17 

16 

1  .00 

45. 

2000. 

17 

18 

17 

0.50 

45. 

2000. 

10 

19 

10 

1  .60 

45. 

2000. 

19 

20 

19 

1  .20 

45. 

2000. 

20 

21 

20 

1  .00 

45. 

2000. 

16 

22 

21 

1.10 

45. 

2000. 

22 

23 

22 

1  .60 

45. 

2000. 

24 

23 

23 

2.70 

45. 

2000. 

24 

25 

24 

2.10 

45. 

2000. 

25 

26 

25 

1  .  10 

45. 

2000. 

26 

27 

26 

1  . 10 

45. 

2000 . 

23 

22 

27 

1  .40 

45. 

2000. 

28 

29 

20 

1  .20 

45. 

2000. 

29 

30 

29 

2.20 

45. 

2000. 

28 

9 

30 

1.10 

45. 

2000. 

9 

31 

31 

1.10 

45 . 

2000. 

31 

32 

32 

1.10 

45. 

2000. 

32 

21 

33 

2.30 

45. 

2000 . 

33 

34 

34 

2.20 

45. 

2000. 

33 

8 

35 

1  .  10 

45. 

2000. 

8 

35 

36 

1.10 

45. 

2000. 

35 

36 

37 

1  . 10 

45. 

2000. 

36 

20 

33 

2.90 

43. 

2000. 

20 

37 

39 

3  80 

45. 

2000. 

19 

38 

40 

3.00 

45. 

2000. 

6 

18 

41 

1  .70 

45. 

2000. 

18 

39 

42 

1.60 

45. 

2000. 

39 

40 

43 

0.90 

45. 

2000  . 

40 

41 

44 

1  .70 

45. 

2000. 

39 

18 

45 

1  .60 

45. 

2000. 

17 

42 

40 

1  . 10 

45. 

2000. 

42 

39 

47 

1.90 

45. 

2000. 

43 

24 

4C 

1  .70 

45. 

2000. 

24 

44 

49 

0.60 

43. 

2000. 

44 

45 

30 

2.50 

43. 

2000. 

43 

46 

51 

0.80 

45. 

2000. 

46 

47 

52 

1  . 10 

45. 

2000. 

40 

49 

53 

1  .60 

45. 

2000. 

49 

50 

34 

1  .20 

45. 

2000. 

50 

5 

55 

3.00 

45. 

150. 

5 

16 

56 

1  .40 

45. 

2000. 

16 

51 

57 

1  .80 

43. 

2000. 

51 

52 

50 

1  .70 

43. 

2000. 

52 

53 

59 

2.10 

45. 

2000. 

52 

40 

60 

3.30 

35. 

2000. 

54 

55 

«1 

2.20 

35. 

2000. 

35 

1 

62 

1 .30 

45. 

2000. 

23 

56 

63 

3.10 

45. 

2000. 

56 

47 

64 

1  .70 

45. 

2000. 

47 

57 

65 

1  .70 

43. 

2000. 

50 

59 

66 

1 .20 

43. 

2000. 

59 

60 

67 

1 .30 

43. 

2000. 

60 

61 

68 

0.70 

45. 

2000. 

58 

62 

69 

0.70 

45. 

2000. 

62 

63 

70 

0.90 

45. 

2000. 

63 

27 
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Table  E-2.  (Continued.) 


LENGTH 

TREE  FLOV 

CAPACITY 

START 

END 

NUMBER  OF 

LINK  I D 

(MILE) 

SPEED (MPH) 

(VEHICLEs/HRl 

Nnnr 

NOnF. 

IJfcNFS 

71 

1 .90 

45. 

2000. 

27 

64 

i 

72 

1 .60 

45. 

2000. 

64 

65 

i 

73 

3.60 

45. 

2000. 

65 

66 

i 

74 

0.40 

45. 

2000. 

67 

68 

i 

75 

1  .00 

45. 

2000. 

68 

15 

i 

76 

1  .40 

45. 

2000. 

15 

59 

i 

77 

1  .00 

45. 

2000 . 

68 

62 

i 

78 

1  .30 

45. 

2000. 

63 

69 

i 

79 

1  .20 

45. 

2000. 

69 

70 

i 

89 

2.20 

45. 

2000 . 

70 

71 

i 

81 

1  .20 

45. 

2000. 

64 

70 

i 

02 

3.90 

45. 

2000. 

25 

65 

i 

83 

2.20 

55. 

2000 . 

72 

73 

i 

84 

1  .00 

55. 

2000. 

73 

61 

i 

85 

2.00 

55. 

200. 

61 

74 

i 

86 

2.60 

55. 

2000. 

72 

75 

i 

87 

1  .70 

55. 

2000. 

75 

76 

i 

88 

1  .40 

55. 

2000 . 

77 

76 

i 

89 

3.30 

55. 

2000. 

77 

78 

i 

90 

1.50 

55. 

2000. 

70 

79 

i 

91 

1  .70 

55. 

2000 . 

1 

80 

i 

92 

1  .40 

55. 

2000. 

80 

81 

i 

93 

1  .70 

55. 

2000. 

81 

82 

i 

94 

2.70 

55. 

2000. 

82 

83 

i 

95 

3.60 

55. 

2000  . 

83 

84 

i 

96 

1  .30 

55 . 

2000. 

84 

79 

i 

97 

2.70 

55. 

2000. 

79 

85 

i 

98 

1  .  io 

55. 

2000. 

85 

86 

i 

99 

0.90 

45. 

2000. 

1 

87 

i 

190 

3.50 

43. 

2000. 

87 

22 

i 

10! 

1  .80 

45. 

2000. 

22 

51 

i 

102 

2.20 

55. 

2000. 

76 

88 

i 

103 

1  .70 

55. 

2000  . 

88 

89 

i 

104 

1  .20 

55. 

2000. 

90 

91 

i 

105 

1  .30 

55. 

2000. 

91 

92 

i 

106 

1  .60 

55. 

2000. 

92 

86 

i 

107 

1  .20 

55. 

2000. 

90 

93 

i 

108 

0.7© 

55. 

2000. 

93 

94 

i 

109 

2.10 

55. 

2000. 

94 

95 

i 

1  10 

2.30 

55. 

2000. 

95 

86 

i 

1  1 1 

3.40 

55. 

2000. 

94 

96 

i 

1  12 

3.20 

55. 

2000. 

96 

97 

i 

1  13 

3.40 

55. 

2000 . 

97 

98 

i 

1 14 

1  .00 

55. 

2000. 

98 

30 

i 

115 

0.80 

45. 

2000. 

99 

34 

i 

1 16 

1  .70 

45. 

2000. 

34 

30 

i 

1 17 

0.70 

45. 

2000. 

30 

100 

i 

1  18 

2.10 

45. 

2000. 

101 

10 

i 

119 

1  .00 

45. 

2000. 

28 

101 

i 

120 

1  . 10 

45. 

2000. 

35 

31 

i 

121 

1.10 

45. 

2000. 

36 

32 

i 

122 

1.10 

45. 

2000. 

102 

103 

i 

123 

0.70 

45. 

2000. 

103 

100 

i 

i 


f 


i 


Table  E-3.  Input  parameters  of  the  sample  problem  for  the  Rancho  Seco  Nuclear  Power  Plant 


INPUT  PARAMETERS 


PEOPLE  PER  VEHICLE  ■ 
VEHICLE  LENCTH(MI)  = 
TIME  I NCREMENT ( HR )  = 
NUMBER  OF  TIME  INCREMENTS  = 
REI'CRTINC  T I  ME  ( HR ) 

NUMBER  OF  NODES  « 
NUMBER  OF  LINKS  * 
CUTOFF  FRACTION  * 


2.0 

4.00E-03 

0.010 

1000 

0.250 

103 

123 

0.9000 


SPECIFIC  PARAMETERS  FOLLOW i 

NOTIFICATION  TIME1HR)  *  0.25 

PREPARATION  TIME(DR)  =  0.00 
DISPERSION  UPDATE! HR )=  0.25 

NUMBER  OF  EVACUATION  ZONES  =  1 


ZONE  NUMBER  1  IS  OUT  TO  10.0  MILES 
INCLUDING  SECTORS: 

1  2  3  4  5  6  ?  8  9  10  11  12  13  14  15  16 
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